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Several ideas for solving the problem of fermion mass hierarchy and mixing and spe- 
cific supersymmetric models that realize it are reviewed. In particular, we discuss many 
models based on SO(10) in four dimensions combined with a family symmetry to ac- 
commodate fermion mass hierarchy and mixing, including the case of neutrinos. These 
models are compared and various tests that can be used to distinguish these models are 
suggested. We also include a discussion of a few SO(10) models in higher space-time 
dimensions. 
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1. Introduction 

The origin of fermion masses and mixing and CP violation is the least understood 
aspect of the Standard Model (SM) of particle physics. In the SM, the fermion 
masses and mixing angles are co mpletely arbitrary. In order to accommodate their 
diverse values (see Tabled) P^EItne Yukawa couplings must range over five orders of 
magnitude. Another relevant aspect is that neutrinos are massless in the framework 
of SM, but recent experiments strongly indicate that they do have small but non- 
vanishing masses. Incorporating these into theory leads to an increase in the number 
of parameters. The flavor problem thus consists the following aspects: how to reduce 
the number of parameters in the Yukawa sector, how to obtain an explanation of 
the mass hierarchy, how to obtain small neutrino masses and large leptonic mixing 
angles. 

As we extend the SM to the Minimal Supersymmetric Standard Model (MSSM), 
the particle spectrum is doubled and many more parameters are introduced into 
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Fig. 1. Diagram of AS = 2 transitions. 



the model. One thus expects large flavor changing neutral current (FCNC) due 
to the presence of squarks. The strongest constraints come from the lighter two 
generations due to the type of diagram given in Fig^ 

^^<6xl0-3(^)(sin^ (1) 

d 

<Z< < w -^^f (2) 

where 9 and 9' refer to the mixing angles in the squark and slept on sectors in the 
bases where the mass matrices of (d s) and (e /d) are diagonal™"' Thus, in order to 
suppress the supersymmetric (SUSY) FCNC, a near degeneracy between the first 
and the second generations of squarks and sleptons is required. This is sometimes 
called SUSY flavor problem. 

There have been many SUSY models proposed to accommodate the observed 
masses and mixing angles. These models can be classified according to the family 
symmetry implemented in the model. We also discuss other mechanisms that have 
been proposed to solve the problem of the fermion mass hierarchy and mixing. 

In the following two sections, we introduce various tools that have been used 
to solve the flavor problem in quark sector and lepton sector, respectively; in Sec. 
El 50(10) is reviewed. This is then followed by a review in Sec. El and Sec. El on 
various models based on 50(10) in AD, combined with mass texture ansatz and 
family symmetry, respectively; a comparison of these models is given at the end of 
Sec. El Sec.[7|is devoted to models based on 50(10) in higher space-time dimensions. 
Sec. El concludes this review. 

Table 1. Current status of fermions masses and CKM matrix elements at M z . 

m u = 1.88 - 2.75 MeV, m c = 616 - 733 MeV, m t = 168 - 194 MeV 
m,/md = 17 - 25, m s = 80.4 — 105 MeV, m b = 2.89 - 3.11 GeV 

m e = 0.487 MeV, m u = 103 MeV, m T = 1.75 GeV 



' 0.9745 - 0.9757 0.219 - 0.224 0.002 - 0.005 \ 
\V C KM,ex P \=\ 0.218 -0.224 0.9736 - 0.9750 0.036-0.046 

0.004 - 0.014 0.034 - 0.046 0.9989 - 0.9993 / 
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2. Quark Masses and Mixing 

Because the strong interaction eigenstates (same as the mass eigenstates) and the 
weak interaction eigenstates do not match, flavor mixing arises. However, currently 
we do not have any fundamental understanding for such a mismatch between the 
strong eigenstates and the weak eigenstates. In the weak basis, the quark mass 
terms are 

Cmass = -Y U U R Q L H U - Y d D R Q L H d + h.c. (3) 

where Q stands for the SU(2) doublet quark; U and D are up- and down- type 
fermion SU (2) singlet; H u and Hd are Higgs fields giving masses to up- and down- 
type quarks. The charged current interaction is given by 

£ cc = A={W+U Llll D L ) + h.c. (4) 

The Yukawa couplings Yi, (i = u,d,e,VLn) are in general non-diagonal. They are 
diagonalized by the bi-unitary transformations 

Yf aB = V UR Y u Vl = diag(y u , y c , y t ) (5) 
Yd™ 9 = v dn Y dV ] dL = diag{y d , y s , y b ) (6) 

where Vr and Vj, are the right-handed and left-handed rotations respectively, and 
all the eigenvalues j/,'s are real and non- negative, and are obtained by diagonalizing 
the hermitian quantity Y^Y and YY^ . The Cabbibo-Kobayashi-Maskawa (CKM) 
matrix is then given by^ 

(V u d v us V u b\ 
Vcd V cs V cb . (7) 
Vd V ts V tb J 

The unitary matrix Vqkm has in general 6 phases. By phase redefinition of various 
quark fields, one can remove 5 of the 6 phases. The remaining one phase is one of the 
sources for CP violation in the quark sector. There are many ways to parameterize 
the CKM matrix, for example, 

r i r i q i r i „9 p -iS q 

C 12 L 13 i '12 c 13 *13 e 

VcKM = | — s 12 c 23 ~ c 12 s 23 s 13 el 9 c 12 c 23 — s 12 s 23 s 13 el 9 s 23 c 13 I ■ (^) 

„i „? _ „i „i „q p iS q _„q q _ Q q „q q p i& q „q q 

, 12''23 L 12 L 23 i 13 e L 12 6 23 *12 L 23''13 e C 23 L 13 

Defining A = s\ 2l AX 2 = and AX 3 (p — iff) = sf 3 e _l<5<J , we obtain an alternative 
parameterization, the Wolfenstein parameterization^ 

/ 1 - \\ 2 X 
Vckm = -A 1 - \X 2 AX 2 | . (9) 

\A\ 3 (1- p-irj) -AX 2 
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Fig. 2. The Unitarity triangle. 

Here the parameters A, p, r\ are of order 1, and A is the sine of the Cabbibo angle 
which is about 0.22, and thus a good choice as an expansion parameter. A param- 
eterization independent measure for the CP violation is the Jarlskog invariant^ 
defined as 

j q CP = im{v ud v; s v; d v cs }. (10) 

Unitarity of Vckm applied to the first and the third columns leads to the following 
condition 3, 

v u v; h + v ud v: b + v cd v: h = Q. (ii) 

A geometrical representation of this equation on the complex plane gives rise to the 
CKM "unitarity triangle" shown in Fig|2 The three angles of the CKM unitarity 
triangle are 

a = Arg{-g#L), ^Aj(-S, ^Arg{-^-). (12) 
VudV ub VtdV tb VcdVc 

The Jarlskog invariant J^ P is proportional to the area of the unitarity triangle. A 
non-vanishing value for J CP thus implies non-vanishing values for (a, j3, 7) which 
in turn indicates CP violation. 



2.1. Textures of Mass Matrices 

There have been many mass textures, with different elements having zeros in the 
mass matrices, proposed in order to accommodate the observed fermion mass hier- 
archy and mixing pattern. Imposing texture ansatz on mass matrices reduces the 
number of parameters in the Yukawa sector; as a consequence, masses and mixing 
angles may be related in some simple ways. This is illustrated in a simplified two 
family example in the up- and down-quark sectors, in which the mass matrices of 
the up- and down-type quarks are assumed to be symmetric and each contains one 
zero entry, 

*-(::)■ M °={ii)- <"» 

a There are five other similar conditions one can write down. 
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Table 2. Mass Texture combinations for up- and down-type quarks with five zeros proposed by 
Ramond, Roberts and Ross. 

I II III IV V 

/ A u \ / A u \ / E u \ / A u \ / E u \ 
M u \ A u B u \ \ A u C u \ B u A u B u C u B u C u ) 

V D u J \ C u D u J \ E u D u J \ Cu DuJ \E U C u D u J 

I A d \ / A d \ / A d \ / A d \ / A d \ 
M d A d B d C d \ A d B d C d \ \ A d B d C d \ \ A d B d A d B d 

V C d D d / V C d D d J V C d D d J V D d J \ D d J 



In this case, there is a very simple relation between the Cabbibo angle and the 
quark masses, 

| Ks | = | /^_ e - S (M) 
V m s V TO c 

where the CP violating phase a arises as the relative phase that enters when com- 
bining the up- and down-quark rotaion matrices. This relation is in good agreement 
with experiment. Note that because m u /m c is quite small compared to md/m s , the 
value of \V US \ is not very sensitive to the complex phase a. One can then generalize 
this to consider the three family case. The ultimate goal of studying texture zero 
ansatz is that it may help us to understand the underlying theory of flavor, if the 
zeros are protected by some family symmetry. In the quark sector, assuming sym- 
metric mass matrices, the total number of texture zeros is at most six, because there 
are six different quark masses. Nevertheless, Ramond, Robert and RosPfound that 
the observed masses and mixing angles cannot be accommodated with six texture 
zeros. They found five combinations of five-zero texture for up- and down-type 
quark mass matrices which give rise to predictions that are consistent with current 
observations for fermion masses and mixing angles. One should note that in the 
context of a grand unified theory (GUT), the texture ansatz is valid only at the 
GUT scale. The vanishing entries in the mass matrices at the GUT scale will be 
filled in by radiactive corrections at lower energy scales. These five solutions are 
summarized in Table [21 

Lop-sided textures have also been considered in model building. In section V, 
we classify various SO(10) models according to whether the mass textures in the 
models are symmetric or lop-sided. Symmetric textures arise if S'O(IO) breaks down 
to the SM group with the left-right symmetric group as the intermediate symmetry, 
while lop-sided textures arise if the intermediate symmetry is SU(5). 

2.2. Froggatt- Nielsen Mechanism and Family Symmetry 

A prototype scenario which prod uces hierarchy in the fermion mass matrices is the 
Froggatt-Nielsen mechanismPSl The idea is that the heaviest matter fields acquire 
their masses through tree level interactions with the Higgs fields while masses of 
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Fig. 3. Schematic diagrams for Froggatt-Nielsen mechanism. Here a and b are the family indices. 
(XiX) are the vector-like Froggatt-Nielsen fields. Figure (a): The tree level diagram generating 
the mass of the third family is given; (b): The mass of the lighter matter fields generated by this 
diagram is ~ 0(( < ®j > ) 2 ); (c): Higher order diagrams generate mass ~ 0(( < ®j > ) n ). 

lighter matter fields are produced by higher dimensional interactions involving, in 
addition to the regular Higgs fields, exotic vector- like pairs of matter fields and the 
so-called flavons (flavor Higgs fields) . Schematic diagrams for these interactions are 
shown in FigEl After integrating out superheavy vector-like matter fields of mass 
M, the mass terms of the light matter fields get suppressed by a factor of ^jf-, 
where < 8 > is the VEVs of the flavons and M is the UV-cutoff of the effective 
theory above which the flavor symmetry is exact. When the family symmetry is 
exact, only the (33) entry is non-zero. When the family symmetry is spontaneously 
broken, the zero entries will be filled in at some order 0(^p). Suppose the family 
symmetry allows only the (23) and (32) elements at order O(^jp), 




SSB 




(15) 



Then a second fermion mass is generated at order 0((^j^) 2 ) after the family 
symmetry is spontaneous broken. The fermion mass hierarchy thus arises. 

To illustrate how the Froggatt-Nielsen mechanism works, suppose there is a 
vector-like pair of matter fields (x®X) with mass M and carrying the same quantum 
numbers as %[)r under the vertical gauge group (e.g. SM or SO(10)), but different 
quantum numbers under the family symmetry. It is therefore possible to have a 
Yukawa coupling yxiphH where H is the SM doublet Higgs if the family symmetry 
permits such a coupling. In addition, there is a gauge singlet 9 which transforms 
non-trivially under the family symmetry. Suppose the coupling y ip^x® is allowed 
by the family symmetry, we then obtain the following seesaw mass matrix, upon H 
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and 6 acquiring VEV's 



Diagonalizing this matrix gives the following mass term for ip 

yy <H ><8> 
~ M ' 

So the suppression factor ^j^- is due to the mixture between the light states and 
the heavy states. This is very simil ar to how the light neutrino masses are generated 
in the seesaw mechanism! ^ ^ 2 [13| 

So what are the possible family symmetries that can be incorporated with the 
Froggatt-Nielsen mechanism? The kinetic terms and gauge interactions of the SM 
have a very large family symmetry, [[/(3)] 5 , where the U(3) factors act on the 
right- and left-handed multiplets of quarks and leptons, respectively. If right-handed 
neutrinos are included, the family symmetry becomes [C/(3)] 6 . Therefore, any family 
symmetry group proposed to incorporate the Froggatt-Nielsen mechanism must 
be contained in [[/(3)] 6 . If all the particles in each family arc unified into one 
single multiplet, as in the case of S'O(IO), the maximal possible family symmetry 
group is reduced to J7(3). The family symmetry can either be global or gauged. 
However, in the case of a global symmetry, there are problems associated with 
the massless Goldstone bosons when the symmetry is broken, and with possibly 
large gravitational quantum corrections. These problems do not arise in the case of 
gauged symmetries. In what follows, we discuss separately the abelian group and 
non-abelian group as a family symmetry. 



2.2.1. Abelian Case 

In compactified string theories, one usually obtains aboundant Abelian symmetries 
below the compactification scale, in additional to the SM gauge group (or a GUT 
gauge group). Suppose that the flavon field 8 and 8 carry +1 and —1 charges under 
U(l), and the ratios ^j^- and have approximately the value of the Cabbibo 
angle 0.22 (or equivalently, the parameter A in the Wolfenstein parameterization). 
By assigning different U{\) charges to different family, one in general, obtains a 
mass matrix of the form 

A™ 11 A" 12 A™ 13 \ 

A™ 21 A" 22 A™ 23 . (18) 

yn 3 i ^"32 ynss J 

The exponent is given by Ql.i + Qhj, where Ql,% {Qn,i) is the U(l) charge 
of the left-handed (right-handed) field of the i-th family. It is usually not difficult 
to find solutions for the charge assignments that are consistent with experiments, 
and the solutions are not unique. An interesting model for both quarks and leptons 
based on anomalous U(1)h was proposed by Irges, Lavigna and Ramond which 
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makes use of the Green- Schw artz anomaly cancellation condition to constrain the 
U(l) charge assignments P^l Due to the Abelian nature, models of this type have 
limited predictive power, because only the powers of the small expansion parameter, 
A, are determined while the relative strengths between different entries are left 
un-determined. As a consequence, every entry in the above mass matrix has an 
un-known 0(1) coefficient associated with it. 

Non-anomalous [ 7(1) n has also been considered in a model constructed by Mira, 
Nardi and RestrepoP^ In this model, anomaly cancellation condition leads to a 
massless up-quark, thus solving the strong CP problem. 

2.2.2. Non-Abelian Case Based on SU(2) 

Models with non-Abclian family symmetry generally have more predictive power 
because the relative strengths between different matrix elements can be determined. 
The original motivation of using non- Abelian group as the family symme try is to 
solve the SUSY flavor problem. SU (2) was proposed by Barbieri et al as a 
family symmetry. It has two attractive features: (i) As we have seen previously, the 
constraints from the SUSY FCNC requires that 

m\-m\ < 10~ 3 to 

to 2 ~ sin</> l 300GeV J 1 ' 

where (m| — m\ 2 ) ~ ni 2 is the average scalar mass squared and tf> is some relevant 
CP phase. SU(2) gives rise to the degeneracies between 1-2 families needed to 
suppress the supersymmetric FCNC in the squark sector; (ii) A multi-step breaking 
of SU (2) gives rise to the observed inter-family hierarchy naturally. Unlike models 
based on the U(l) family symmetry, in which one has the freedom in choosing U(l) 
charges for various matter fields, a SU (2) family symmetry appears to be a much 
more constrained framework for constructing realistic models. 

The heaviness of the top-quark and suppression of the SUSY FCNC together 
suggest that the three families of matter fields transform under a SU(2) family 
symmetry as 

ie^ = 2ei (20) 

where a = 1, 2 and the subscripts refer to family indices. In the symmetric limit, only 
the third family of matter fields have non- vanishing masses. This can be understood 
easily since the third family of matter fields have much higher masses compared to 
the other two families of matter fields. SU (2) breaks down in two steps: 

57/(2) u{l) ^ nothing (21) 

with e' <C e <C 1 and M is the UV cut-off of the effective theory mentioned before. 
These small parameters e and e 1 are the ratios of the vacuum expectation values of 
the flavon fields to the cut-off scale. Note that because 



V'aV'a ~ Is, Ma ~ 2, tp a ip b ~2<g)2 = lA83 



(22) 
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the only relevant flavon fields are in the 1a, 2 and 3 dimensional representations of 
SU(2), namely, 



A ab 



s ab 



(23) 



So a generic mass matrix constrained by SU (2) family symmetry is of the following 
form: 



< s >, 


< 1a > 


< <t> > 


< <? 




1 



(24) 



7 



To see how the vacuum alignment in the flavon sector is achieved, let us first 
consider the supersymmetric limit with only one conjugate pair of doublets (0 ©</>), 
anti-symmetric singlets (A © A) and triplets (S® S) of SU(2). The most general 
renormalizable superpotential is then given by^J 



WflavonM 0, S, S) = cbS<j> + 4>S(f> + X^U + X S SS + X A AA (25) 

where X s , X^ and Xa are dimensionful parameters in the superpotential. Note that 
the anti-symmetric singlet fields A and A are decoupled from other fields. From the 
F-flat conditions, one obtains the following solutions, 



E 



< s ab > = - < 4> a >< 4>b > /x s 



(26) 
(27) 



X A A ab = 0. (28) 

Thus the relative strengths of < S > and < (j> > are determined. For Xa ^ 0, the 
F-flat conditions imply < A >=< A >~ 0. Non- vanishing < A > and < A > can 
be obtained if non-renormalizable operators arc introduced!^ 

If all the 16 observed matter fields in one family form a single representation as in 
the case of S'O(IO), the most general effective superpotential, after integrating out 
all the heavy Froggatt-Nielsen fields, that generates fermion masses for a 5*0(10) x 
SU(2) model has the following very simple form 



W 



Cjab 



A ab 



In a specific SU(2) basis, 

M . ^f^\ ( Aab ) 



M 



O 



M 



O 







(s ab ) 

M 



O 



e e e 



(29) 



(30) 



Here we have indicated the VEVs all the flavon fields could acquire for symmetry 
breaking in Ea. (|21|l . The mass matrix M takes the following form 



M ~ O | -e' e e 
e'e e 1 

The hierarchy is thus built into this mass matrix. 



(31) 
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2.2.3. Non-Abelian Case Based on SU(3) 

Ultimately, the SU(2) family symmetry can be embedded into SU(3), under which 
the three families form a triplet. A model based on SU(3) is presented in Ref. 
18. The three families form a SU(3) triplet before the symmetry is broken. The 
symmetry breaking takes place at two steps 

SU(3) — ► SU{2) — ► nothing. 

The SU{3) anti-triplet flavon fields, 3 and 23 , acquire VEV's along the following 
directions, breaking the SU (3) symmetry 




< 3 >= , < 23 >= lb] (32) 



where < 3 > triggers the first stage of breaking, and < 23 > triggers the second 
stage of breaking. The leading order Yukawa couplings of the matter fields are 

H(jL^ 3 ^i + _1_^; 3 ^ 3 ). (33) 



This gives rise to a mass matrix of the following form 

(34) 



/0 

b 2 b 2 



MS „M: 







His 









\ Xr- — — — v- —I— — — 

Assuming M 3 ~ a 3 ^S> M23 3> b and e = 6/^/23, one thus obtains a hierarchical 
structure in the (2,3) block. There are also operators that mix 3 with (j> 23 , 

e 2 — i —i—- 
jj— — ^(^023^03 + ^03^-023)- (35) 

These operators give rise to contributions to the (23) and (32) matrix elements of 
order 0(e), which is larger than the (22) element of order C(e 2 ), leading to a wrong 
prediction for V c t- One way to suppress these operators is to impose a discrete Z2 
symmetry, under which <f> 3 and 4> 23 have opposite parity. Thus the operators given 
in Ea. H33J) are allowed by the Z 2 symmetry while the operators given in Eq.(j2HJ) 
only arise at higher order with a suppression factor e 2 . To fill in the first row and 
column, one has to introduce additional SU(3) triplet flavon fields, <p 3 and 023, 
which acquire VEV's along the following directions 

< 03 >= ( j , < 023 >= \^ bj (36) 

and consider the following operators 

T^(e ijk ^ihzM) H (37) 

lg 4^(e^V^3 J 023,fc)(6 im "V'r03, m 023,„)-ff. (38) 
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The operators given in Ea. (|37|l generate the entries (12) = (13) = —(21) = —(31) 
of order 0(e 3 ). These entries are anti-symmetric due to the presence of the anti- 
symmetric tensor t % ^ k in the couplings. The operators given in Eg. 1(38(1 generate the 
(11) matrix element of order C(e 8 ). The suppression in this operator is due to the 
presence of the Z2 symmetry discussed above and an additional i?-symmetry; these 
symmetries also forbid all the operators which could lead to un-realistic predictions. 
The operators given in Ea. 1(33(1 . 1(35(1 . ((37(1 and 1)38(1 together give rise to a Yukawa 
matrix of the form, in the leading order of the expansion parameter, e, 

'0(e 8 ) Ae 3 Ae 3 \ 
-Ae 3 e 2 e 2 (39) 
-Ae 3 e 2 1 / 

which can accommodate realistic fermion masses and mixing angles. The vacuum 
alignment leading to Eo . I(32|l and 1(36(1 is discussed in detail in Ref.18. 

We comment that the absolute mass scale of the family symmetry, M, cannot be 
determined in the Froggatt-Nielsen type of scenario, because it is the ratio, (e, e ), 
rather than the absolute mass scale, M, that is phenomenologically relevant. Some 
attempts have been made by having the SUSY breaking messenger fields play also 
the role of Froggatt-Nielsen fields such that the family symmetry scale is linked to 
the SUSY breaking scale. Though attractive, these models have difficulties getting 
low SUSY breaking scale. ^ 

2.3. Ideas from Extra Dimensions 

2.3.1. Split Fermion Scenario in Factorizable Geometry 

It has been proposed that the existence of extra compact spatial dimensions could 
account for the large hierarchy betwee n th e Planck scale and the electroweak scale. 
Arkani-Hamed, Dimopoulos and DvalP^l showed that extra dimensions of size ~ 
1 /TeV may provides a solution to the gauge hierarchy probl em. 

Based on this framework, Arkani-Hamed and SchmaltJ^ proposed a non- 
supersymmetric model in which the mass hierarchy is generated by localizing zero 
modes of the weak doublet and singlet fermions at different locations. Note that 
this mechanism works despite of the size of the extra dimensions. Consider a chiral 
fermion iff in 5D. The action for W coupled to a background scalar field <J> is given 

by 

S = J d 4 dy^ [i^d„ + ij 5 d 5 + $(y)] (40) 
The chiral fermion ^>(x,y) can be expanded in the product basis 

V(x,y) =J2< y\ L > n> p L^n{x)+ < y\R,n> P R <if> n (x). (41) 

n 

Here \L,n > and \R,n > satisfy the following equations 

aa^\L,n > = {-dj + $ 2 + $)|i, n >= ^ \L, n > (42) 
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a?a\R, n > = {-d\ + $ 2 - $)\R, n >= (fi\R, n >, (43) 
where $ = 3$$ and 

a = d 5 + <S>(y), t=-ft + *( tf ). (44) 

For a special choice of linearized background field, &(y) = 2^i 2 y where \i is a 
constant of mass dimension one, the operator a and a* become the usual creation 
and annihilation operators of a simple harmonic oscillator. In what follows, we will 
concentrate on this special case. Expanding the 5D action in terms of \L, n > and 
\R,n >, and then integrating out y, we obtain the 4D effective action 

S = / d 4 $ L i^d^P L ij L + ^ R irduP R ipR + !>Jyfdy. + M*MJ, (45) 

^ n=l 

where the first two terms are kinetic terms of the chiral zero modes whose wave 
functions are Gaussian centered at y = 0, 

In general, there is a bulk mass term, to'S'I/, in the 5D action, Ea. (|40() . because mass 
terms for 5-dimensional field are allowed by all symmetries. In this case, instead 
of being centered at y = 0, the Gaussian wave function is centered at y = m/2/i. 
Thus by having different bulk mass term, to, different bulk fields can be localized 
at different AD slices in the 51? bulk. 

For simplicity, the wave function of the Higgs doublet is assumed to have con- 
stant spread along the fifth dimension. In this case, the Yukawa coupling is 



S y uk = vu J d 4 xH(x)f i (x)f^(x) J dy^M^fM, (47) 

where (y) and (j>z c (y) are the zero mode wave functions of the SU(2) doublet, /j, 

and singlet, /J, respectively. Integrating out the y coordinate, the effective Yukawa 
coupling is then given by 

Vii J dy<j> f M<i> f y(y) = Wi^f / dye-^-^V^^ 2 = e~^H\ (48) 

where Tf i = m^/2/j, and rf c = m;c/2/j are the locations at which the Gaussian 
wave functions <j)t i and (pfaiy) are centered, and = |r? ? — TtA is the "distance" 
between the peaks of these two wave functions. Thus the large hierarchy among 
fcrmion masses can be generated by having different TOj. 

Variations based on this mechanism have been investigated in Ref. 22, 23, 24. 
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2.3.2. Non-Factorizable Geometry 

Randall and SundrurrP^l proposed an alternative based on non-factorizable geome- 
try from which the gauge hierarchy can be derived. In their original proposal, only 
gravity propagates in the bulk. If the SM Higgs doublet is confined to the TeV 
brane, which is required if the gauge hierarchy is assumed to arise from the warped 
geometry, while all other SM particles are allowed to propagate in the bulk, it is 
possible to understand the fermion mass hierarchy in this setup. The equations of 
motion for various bulk fields are given in the following compact forirk 

(e^^/A + e s °d 5 (e- s °d 5 ) - M|)$(^, y) = 0, 



26 27 



(49) 



where for <f> 



-2a 



* L , fl ) we have M| = (ak 2 + ba" (y), C(C ± l)fc 2 ± Co" (y)) 
and s — (4, 1), where a and C are bulk mass terms of the scalar and fermionic fields, 
and 6 is a boundary mass term for the scalar field. The field y) is decomposed 

into an infinite sum of Kaluza-Kline (KK) modes as follows, 

1 



$(^y) = -=]>>(„) (^)/„(y). 



The profile of the n-th mode, f n (y), satisfies 
(_ e - 95 ( e — d 5 )+Ml)f n (y) 



„2(T 



(50) 



(51) 



where M% — (ak 2 , C(C ± l)fc 2 )- m n is the mass of the n-th KK mode. The solution 
for the zero modes of the bulk spin-1/2 fields are found to bel22IZII 



My) = w e- ca , 

iV 

where c = C for left-handed fermions and c = — 
normalization constant 1/Nq is given by 

1 (1 - 2c)nkR 



N 2 



e (l-2c)irkR _ I 

Thus the bulk fermion can be decomposed into 



(52) 

C for right-handed fermions. The 

(53) 



9 L (x»,y) = e 2 °<l>(x^y) 



1 



^e (2 - c) ^(0)(^) + --' 



The Yukawa coupling for the charged fermions to the Higgs doublet reads 



Y 

~M~n 



d 4 x 



-gV R ,i(x, y)*L,j(x, y)H(x)S(y - irR), 



(54) 



(55) 



where Y^ are dimensionless 0(1) coefficients. The effective Yukawa coupling in four 
dimensions is obtained after integrating out the fifth coordinate, y: 



M p i y 2ttR 



(1 - 2c RA )irkR 



o(l-2c R i )TTkR _ I 



1/2 



(1 - 2c Ljj )irkR 

g(l— 2cLj)-xkR _ ^ 



1/2 



xR 



dy\ 



-g e^ CR '' )a e ( - 2 - CL ^ a e a 5(y - ttR) 



-tvR 
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= Y v k y/(l-2cfl,i)(l-2c L ,j) e (l^ t -c tJ )TM (56) 

2 M p ; yj ' ( e {l-2c RA )iskR Z 1 )( e (l-2c i , I ) 7r fci? _ ' V j 

Thus by choosing Ci,, and c^i all of 0(1), we can reproduce the observed mass 
hierarchy and mixing. We note that this mechanism is not predictive in the sense 
that it does not reduce the number of parameters; the virtue of this mechanism is 
that the large hierarchy observed in fermion masses arises from parameters all of 
0(1). A configuration that reproduce the observed mass hierarchy has been found 
by Huber and ShafiPl 

3. Lepton Masses and Mixing 

If neutrinos are massive, the mixing arises in the leptonic charged current interac- 
tion, 

C cc = -j=Ul M (W+V L ^E L ) + h.c. (57) 

The leptonic mixing (LM) matrix ESEDB 1 1 ^ 21 Ulm^ is obtained by diagonalizing 
the Yukawa matrix of the charged leptons and the effective neutrino mass matrix, 
assuming neutrinos are Majorana particles, 

ydiag = VeR Y e vl = diag(y e , y„ y T ) (58) 

Mil a L 9 = Vu LL M%[V* L = diag{m Vl ,m V2 ,m U3 ), (59) 

where V ULL is an orthogonal matrix, and it can be parameterized as a product of 
a CKM-like mixing matrix, which has three mixing angles and one CP violating 
phase, with a diagonal phase matrix, 

U L M = V eL Vj LL 

c 12 c 13 s 12 c 13 s 13 e iSl \ /l 

J A — A e l q 1 p iS i A A — q 1 q 1 q 1 p iS ' s l A I I p i ^i L 

A 12 c 23 c 12 A 23*13 e L 12 c 23 *12 6 23*13 e *23 c 13 II 

— A A s l —A q 1 — s l A q 1 p^i A A I \ P i ~¥" 

'12*23 L 12 c 23*13 e c 12*23 *12 c 23*13 e C 23 L 13 / \ e 

(60) 

which relates the neutrino mass eigen states to the flavor eigenstates by 

We > = U evi \U! > +U eV2 \U 2 > +U eV3 \v 3 > (61) 

K > = u i*»i K > +Upu*\v2 > +U^ V3 \v 3 > (62) 
\v T > = U rVl \vi > +U Tl/2 \is 2 > +U TU3 \i>3 > . (63) 

b The LM matrix is sometimes called Pontecorvo-Maki-Nakagawa-Sakata (PMNS) or Maki- 
Nakagawa-Sakata (MNS) matrix. It wa s first discussed, in a two flavor case, by PontecorvJ221 
and by Maki, Nakagawa and Sakat aEOl Th e mixing matrix with 3 flavors was first discussed by 
Lee, Pakvasa, Shrock and Sugawaral^3S21 
c If neutrinos are Dirac particles, Eg. 1591 becomes 

M diag = M ef f y\ = dia g( mui _ im „ 3 ). 
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Note that the Majorana condition, 

C{vj) T = u h (64) 

where C is the charge conjugate operator, forbids the rephasing of the Majorana 
fields. Therefore, we can only remove 3 of the 6 phases present in the unitary 
matrix Ulm by redefining the charged lepton fields and are left with t hree CP 
violating phases in the leptonic sector, if neutrinos are Majorana particles j ^^ l ^4 | 35 ] 
Thus Ulm can be parameterized as a product of a unitary matrix, analogous to 
the CKM matrix which has one phase (the so-called universal phase), Si, and a 
diagonal phase matrix which contains two phases (the so-called Majorana phases), 
«2i and a3i. The leptonic analog of the Jarlskog invariant, which measures the CP 
violation due to the universal phase, is given by 

J l CP = Im{U^U eV3 U; v3 U* eV2 }. (65) 

For the Majorana phases, the rephasing invariant CP violation measures are 

Sx=Im{U eul U* eV3 }, S 2 = Im{U eU2 U* eU3 }. (66) 

From S\ and S*2, one can then determine the Majorana phases 

cosa3i=1 ~ 2 w^mj 2 (67) 

co S (a 31 -a 21 ) = l-2^-^- F . (68) 

The recently reported measurements from KamLAND reactor experiment!^! 
confirm the large mixing angle (LMA) soluti on to be the unique oscillation solution 
to the solar neutrino problem at 4.7 cr The global analysis including 

Solar + KamLAND + CHOO^S! 

indicate the following allowed region at 3<rpHI 

5.1 x 10~ 5 < Amfi < 9.7 x lO^eF 2 (69) 
0.29 < tan 2 12 < 0.86. (70) 

The allowed regions at 3<r level based on a global fit including SK-*^ 1 + Solar + 
CHOOZ for the atmospheric parameters and the CHOOZ angle ar 

1.4 x 10~ 3 < Ato 2 2 < 6.0 x 10~ 3 eV 2 (71) 
0.4 < tan 2 9 23 < 3.0 (72) 
sin 2 13 < 0.06. (73) 

And the magnitudes of Ulm elements at la (3er) are given by-^ 

/ (0.73)0.79 - 0.86(0.88) (0.47)0.50 - 0.61(0.67) - 0.16(0.23) 
\U L m\ = (0.17)0.24-0.52(0.57) (0.37)0.44-0.69(0.73) (0.56)0.63-0.79(0.84) 
\ (0.20)0.26 - 0.52(0.58) (0.40)0.47 - 0.71(0.75) (0.54)0.60 - 0.77(0.82) j 
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To see what a bi-large mixing means, let us assume, to a good approximation, 
that #i3 = r) <C 1. In this case, the LM matrix reads, 

/ C12 S12 i] \ 

Ulm = 75 vl i/v ^ i ( 75 ) 

(£12— gClg) (C12+??S12) / 

\ -\/2 \/2 / 

which in turn implies the neutrino mass eigenstates given in terms of flavor eigen- 
states as 

1 r) 
\vi > = c 12 \v e > -—j=si 2 { > - \v T > ) - —^ c ^( \ V H > + Wr > ) (76) 

1 Tj 
\V2 > = SuWe > +-y=Cl2{ |^ > - \V r > ) - —f^ S 12( > + I V r > ) (77) 

b > = ^={W„> + Wr > ) + vWe > ■ (78) 

3.1. Generation of Small Neutrino Masses 
3.1.1. Small Neutrino Masses from See-saw Mechanism 

The observation of neutrino oscillations provides the first indication of beyond the 
Standard Model physics. It has two implications: neutrinos have non-zero masses, 
and lepton family numbers are violated. In the SM, neutrinos are massless because 
there are no SU (2) singlet neutrinos nor are there SU (2) triplet Higgss. Adding one 
SU (2) singlet neutrino for each family is the simplest way to introduce neutrino 
masses. Because the right-handed neutrinos are SM singlets, the Majorana mass 
terms for the SU (2) singlet neutrinos are not forbidden by the symmetry. In the 
Lagrangian, there are Dirac m ass term for the neutrinos and the right-handed 
Majorana mass terms IHH2II31 

1 _ T _ 

^seesaw = ~M VLR V R V L - -M VRR V R V R + h.C. 

(In general, there could be a non-vanishing mass term for (ulvl)', this is the Type 
II see-saw mechanismP2l Such a mass term can be obtained in 5*0(10) from the 
Yukawa coupling given in Ea. l|175|l ). After integrating out the heavy right-handed 
neutrinos, 

dC _ T 

- t=- = M LR v L + M RR v R 
ovr 

v R = ^ T Ml R M R l (80) 
one then obtains the effective light neutrino Majorana mass matrix 

M ll = M t Vlr M- rr M Vlr . (81) 
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For Mlr of the order of the weak scale as the mass scale of other charged fermions, 
the Majorana mass term Mrr must be around 10 12 ~ 14 GeV to give rise to neutrino 
masses of the order of 0.1 eV. As we will see later, most grand unified theories 
predict the existence of the right-handed neutrinos. Furthermore, the GUT scale 
provides an understanding why Mrr is large. 



3.1.2. Small Neutrino Masses from Large Extra Dimension 

An interesting way to generate small Dira c neut rino masses arises in models with 
large extra dimensions of size ~ 

(l/TeV).E5EEI Consider the case in which only 
gravity can propagate in the bulk, while the SM particles and interactions are 
confined to the brane. Because the right-handed neutrinos are SM singlets, they 
are the only particles that can propagate in the bulk. Its Yukawa coupling to the 
charged lepton L(x), and the Higgs doublet, H(x), which are confined to the brane, 
is given by 



S = yJd 4 xJ dyL{x)H{x)v R (x, y)5(y). (82) 
Compactificd on a circle 5 1 , vji(x^y) can be decomposed into 

1 oo 

v R {x, y) = v ( R\xV ny/R - (83) 

Below the compactification scale, we thus obtain a Dirac neutrino mass 

y<H> y<H>M™ 



2wRM^ D M i 



4D 



(84) 



where tha last steps follows from the relation between the AD Planck scale and 5D 
Planck scale, 2irRM%f> = ( M 4 f /M^f ) 2 . With M^ D - (1 - 10)TeV which could 
avoid the gauge hierarchy problem, one obtains a highly suppressed m„ which is 
consistent with the experimental observations. 



3.1.3. Small Neutrino Masses from Warped Geometry 

Small neutrino masses of the Dirac type are possible if right-handed neutrinos are 
localized toward the Planck brane while the lepton doublets are localized toward 
the TeV brane. This results in a small overlap between the zero mode profiles of 
the lepton doublets and the right-handed neutrinos based on the formulation given 
in Sec. II. Models of this type have been constructed in Ref. 26, 28, 43. 
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3.2. Bi-Large Neutrino Mixing Angles 



To obtain the bi-large mixing pattern for the neutrinos, in additional to having the 
hierarchical mass pattern, let us consider, for example, the following mass textureP^ 




(85) 

with t < 1 which is a special case of the following texture 

(86) 




first proposed in Ref. 49 in which the elements in (2, 3) block are taken to have equal 
strengths to accommodate near bi-maximal mixing. The modification of adding the 
term t n in the (23) and (32) entries is needed in order to accommodate a large, but 
non-maximal solar angle in the so-called "light side" region (0 < 9 < 7r/4)P"3ft j s 
possible to obtain the LMA solution at 3cr level with n ranging from 1 to 2. 

An interesting alternative in which a 3 x 2 neutrino Dira c mass matrix is consid- 
ered was proposed recently by Kuchimanchi and Mohapatr A 3 x 2 neutrino 
Dirac mass matrix arises if there are only two right-handed neurtinos, instead of 
three. The existence of two right-handed neutrinos is required by the cancellation 
of Witten anomaly, if a global leptonic SU (2) family symmetry is imposecEES 
Along this line, Frampton, Glashow and Yanagida proposed a model, which has the 
following Lagrangianj^l 

C = l( N ^( A M 2 ) (nI) +{NiN2) I I I \ H + h - c - ( 8 7) 

The effective neutrino mass matrix due to this Lagrangian is obtained, using the 
see-saw formula, 




/ a 2 aa n \ 

; Mi u 1 

aa a 2 i b 2 bb 

Mi Mi ' M 2 M s 

V o ^- ^- I 

\ u M 2 M 2 / 



(88) 



where a, b, b are real and a = \a \e lS . By takinging all of them to be real, with the 
choice a — V2a and b — b , and assuming a 2 /Mi <C b 2 /M2, the effective neutrino 
masses and mixing matrix are obtained 

2a 2 2b 2 
m Vl - 0, m V2 = — , 771^3 = — (89) 

'1/V21/V2 \ /l 

U= ( -1/2 1/2 1/V2 x cosfl sin6> ) , (90) 
1/2 -1/2 1/V2/ \0-sm^cos( 
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where 9 ~ m U2 / \f2m V3 , and the observed bi-large mixing angles and Am^ tm and 
ArriQ can be accommodated. An interesting feature of this model is that the sign 
of the baryon number asymmetry (B oc £b = Y 2 a 2 b 2 sin 25) is related to the sign 
of the CP violation in neutrino oscillation (( OS(: ) in the following way 

f- = -]j^( 2 + ya )&«- fl ( 91 ) 

assuming the baryon number asymmetry is resulting from leptogenesis due to 
the decay of the lighter one of the two heavy neutrinos, Ni, whose mass is of 
O(10 10 GeV). A 50(10) model which gives rise to the neutrino mass ansatz, 
Eq. (|87|l , has been constructed!^ A more detailed discussion on this model is given 
in Sec. 6. 

Other phenomenologically viable textures for neutrino mass matrix are analyzed 
in Ref. 55. 



3.2.1. 5O(10) GUT realization 

In 50(10) models in AD, the bi-large mixing in the leptonic sector arises in two 
ways (A detailed classificatio n a ccording to how the maximal — v T mixing arises 
is given by Barr and DorsneJ^j): 

(i) Symmetric mass textures for the charged fermions: 

This scenario is realized in symmetric textures arising from left-right symmetric 
breaking chain of 50(10). In this case, both the large solar mixing angle and the 
maximal atmospheric mixing angle come from the diagonalization of the effective 
neutrino mass matrix. A characteristic of this class of models is that the predicted 
value for | U el/3 | element tends to be larger than the value predicted by models in 
class (ii) below. 

(ii) Lop-sided mass textures for charged fermions: 

In this scenario, the large atmospheric mixing angle comes from charged lepton 
mixing matrix. This scenario is realized in models with SU (5) as the intermediate 
symmetry which gives rise to the so-called "lop-sided" mass textures, due to the 
SU(5) relation, 

M e = Mj. (92) 

Due to the lop-sided nature of M e and Mj, the large atmospheric neutrino mixing is 
related to the large mixing in the (23) sector of the RH charged lepton diagonaliza- 
tion matrix, instead of V c b- It thus provides an explanation why the small value of 
V c b and the large value of U^ 3 exist simultaneously. The large solar mixing angles 
comes from the diagonalization matrix for the neutrino mass matrix. Because the 
two large mixing angles come from different sources, the constraint on U eV3 is not as 
strong as in class (i). In fact, the prediction for f7 e „ 3 in this class of models tend to 
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be quite small. On the other hand, this mechanism also predicts an enhanced decay 
rate for the flavor-violating process, fj, — > e 7, which is close to current experimental 
limit. 

We will discuss these two classes of models in detail in SecEl 
3.2.2. Large Mixing from Renormalization Group Evolution 

It has been shown that it is possible to obtain large neutrino mixing angles through 
the renormalization group Recently, Mohapatra, Parida 

and Rajasekaran observed in Ref. 62 that bi-large mixing angles can be driven by 
the renormalization group evolution, assuming that the CKM matrix and the LM 
matrix are identical at the GUT scale, which is a natural consequence of quark- 
lepton unification. The only requirement for this mechanism to work is that the 
masses of the three neutrinos are nearly degenerate of the form m 3 > m 2 > mi 
and have same CP parity. The one-loop RGE of the effective left-handed Majorana 
neutrino mass operator is given b y I I I I 

—T£- = -{n u m u + m v P + P T m u }, (93) 

where t = ln/x. In the MSSM, P and k u are given by, 

1 YjY e 1 hi 

32tt 2 cos 2 ~ 32?r 2 cos 2 f3" 

T^lU + 6 52 2 ~ 6^f^] ~ -^- 2 [\ 9 \ + 6 52 2 - 6-%], (<■:,) 
I67H 5 sin (3 I67H 5 sm (3 



P = -^2^ ^ -^2^b^(0,0,l) ^ diag(0AP T ) (94) 



where g\ = |gy is the t/(l) gauge coupling constant, Y u and Y e are the 3x3 Yukawa 
coupling matrices for the up-quarks and charged leptons respectively, and ht and 
h T are the t- and r- Yukawa couplings. One can then follow the "diagonal-and-run" 
procedure, and obtain the RGE's at scales between Mr > fi> MgusY for the mass 
eigenvalues and the three mixing angles, assuming CP violating phases vanish, 

= -4P T miU? v . - m t n u , (i = 1, 2, 3) (96) 



dt 
d S23 
dt 

d S13 

dt 
d su 
dt 



= -2P T c 2 2Z {-s 12 U TVl V31 + Ci 2 U TV2 V 32 ) (97) 
= -2P T c 23 cl 3 (c 12 U TUl V 31 + s 12 U TU2 V 32 ) (98) 
= -2P T c 12 (c 23 s 13 s 12 U TUl V '31 - c 23 si 3 ci 2 U TV2 V 32 + U TVl U TV2 V 2 i), (99) 



where Vy = (m t + mj)/ {mi — mj). Because the LM matrix is identical to the CKM 
matrix, we have, at the GUT scale, the following initial conditions, s\ 2 ~ A, s 23 ~ 
(D{\ 2 ) and Si 3 ~ C(A 3 ), where A is the Wolfenstein parameter. (Note that the RG 



d Note that some of the earlier results were not entirely correct; re-derivation of these results has 
been done in Ref. 65, 66. 
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evolution has negligible effect on the Wolfenstein parameter, see Ea. ll38| l). When 
the masses m, and uij are nearly degenerate, Vy approaches infinity. Thus it drives 
the mixing angles to become large. Starting with the values of {m\,m\,w^) = 
(0.2983,0.2997,0.3383) eV at the GUT scale, the solutions at the weak scale for 
the masses are (mi,m 2 ,m 3 ) = (0.2410,0.2411,0.2435) eV, which correspond to 
Am 2 tm = 1.1 x 10~ 3 eV 2 and Atoq = 4.8 x 10~ 5 eV 2 . The mixing angles predicted 
at the weak scale are sin #23 = 0.68, sin #12 = 0.568 and sin 6*13 = 0.08. Because 
the masses are larger than 0.1 eV, they are testable at the present searches for the 
neutrinoless double beta decay. 

Models based on GUT and horizontal symmetry often suffer from fine-tunning 
or the difficulty of constructing a viable scalar potential that gives rise to the re- 
quired vacua. Along the line discussed in the above paragraph, some attempts have 
been made to s how that the maximal m ixing anglt^ ^ Sal anc | nearly degenerate 
neutrino n^WmW4Wiffl are manifestations of infrared fixed points 
(IRFP) of the RGEs given above in Ea. l|96(l - (|99[1 . under certain assumptions. 



3.2.3. Bi-large Mixing and b — r Unification 

In the minimal £0(10) model utilizing Type II see-saw mechanism with one 10 and 
one 126, we have the following relations: 



M u = f < 10 > +h< T26 + > (100) 

M d = / < 10 > +h < 127T > (101) 

M e = / < 10 > -3/i < 126" > (102) 

m »lr = / < 10 > -3h < T26 + > (103) 

where / and h are Yukawa matrices; the mass terms M Vj ll and M Vj rr are both 
due to the coupling to 126, 

M V>LL =h< I26' + > (104) 

M^rr =h< T26'° > (105) 

where the superscripts +/ — /0 refer to the sign of the hypercharge Y (see Table 3). 
The small neutrino masses are explained by the Type II see-saw mechanism with the 
assumption that the LH Majorana mass term dominates over the usual Type I see- 
saw term, thus it is proportional to the Yukawa matrix h, which can be determined 
by calculating the difference between Md and M e . Using down- type quark masses, 
charged lepton masses, and CKM matrix elements, which have roughly the form 

/A 3 A 3 A 3 \ 

M bT - A 3 A 2 A 2 m br , (106) 
\A 3 A 2 1 / 

predictions for neutrino masses and LM matrix elements have been made.^l^ 
The large atmospheric mixing results from a small deviation of C(A 2 ) from b — r 
unification. 
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Predictions from a detailed numerical study made in Ref. 79 are sin 2 26*23 < 0-9 
and sin 2 2#i2 > 0.9, which are experimentally allowed only at 3a level; these unique 
predictions can thus be used to test this type of models. Note that the best fit 
value of sin 2 2#23 cannot be accommodated in these models. The prediction for U e 3 
is about 0.16, very close to the sensitivity of current experiments. We also note that 
as this type of models do not address the origin of the flavor structure, they are 
not as predictive as 5*0(10) models combined with family symmetry, in which as 
non- minimal Higgs content is usually present (see Sec. 6). 

3.3. CP violation in Leptonic Sector 

As we mentioned previously, if neutrinos are Dirac particles, there is only one phase 
in the LM matrix. On the other hand, if neutrinos are Majorana particles, which is 
the case if see-saw mechanism is implemented to give small neutrino masses, there 
are two additional phases. These two types of phases have very different impacts 
on phenomenology. The universal phase, Si , affects the transition probability in the 
neutrino oscillation 

P{v a -> vfi = <W -Aj2Re{U m U fjj U* aj U^} sh 2 (Am^) 

i>j 

+2 ]T ImiUciUuU^Ufc} sin 2 (Am 2 A). (107) 

i>3 

Note that the only chance that one might observe CP violation in neutrino oscil- 
lation is to have LMA solution in the solar sector, and to have large value for #13. 
The Majorana phases affect the matrix element for the neutrinoless double beta 
(PPov) decay, | < m > |, given in terms of the rephasing invariant quantities by 

I < m > | 2 - m\\U el \ A + m 2 |[/ e2 | 4 + m 2 |[/ e3 | 4 
+2m 1 m 2 \U el \ 2 \U e 2\ 2 cos a 2 i 
+2mim 3 |f/ e i| 2 |C/ e3 | 2 cosa 3 i 

+2m2m 3 |t/ e2 | 2 |[/ e 3| 2 cos(a 3 i - a 21 ). (108) 

The current bound on | < m > I from Heidelberg-Moscow experiment is 0.11 — 
0.56 eV at 95% confidence levelPB 

4. Grand Unified Theories Based on SO(10) 

The smallest GUT group SU (5) in its minimal form is very strongly constrained due 
to the non-observation of proton decay!^S2IHSI The next candidate is the rank-5 
5*0(10), which is a very attractive candidate as a GUT group for many reasons: First 
of all, all of its irreducible representations are free of anomaly, unlike in the case of 
SU(5) where the representations of the matter fields, 5 © 10, are carefully chosen 
to cancel the anomaly. It unifies all the 15 known fermions with the right-handed 
neutrino for each family into one 16-dimensional spinor representation. The seesaw 
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mechanism then arises very naturally, and the small but non-zero neutrino masses 
can thus be explained, as evidenced by recent atmospheric neutrino oscillation 
data from Super-Kamiokande indicating small non-zero neutrino masses. Because 
a complete quark-lepton symmetry is achieved, it has the promise for explaining 
the pattern of fermion masses and mixing. In some 50(10) models, R-parity is 
conserved automatically at all energy scales. This is to be contrasted with MSSM 
and SUSY SU (5) where R-parity must be imposed by hand. Because (B — L) is a 
gauge symmetry contained in 50(10), it has the promises of baryogenesis. In what 
follows, we briefly review the structure of 50(10) models. A detail discussion can 
be found in Ref. 84. 

4.1. The Algebra of SO (2n) 

It is convenient to discuss the 50(2ro) algebra in the SU(n) basis P^l Consider a set 
of n operators (i — 1, ...,n), and their hermitian conjugates, Q, satisfying 

{&,£]}=<%, =0, (109) 

where { , } denotes an anti-commutator and [ , ] denotes a commutator. The 
operators K!j defined as 

K) = (HO) 
satisfy the algebra of the U (n) group 

[K i j ,K™]=8™K i n -8i l K™. (Ill) 
We can then define the following 2n operators, T M (fi= 1, 2n) 

r 2j -i = -i ( & - 4 ) 

r 2 , = + j = l,...,n. (112) 

The T M form the Clifford algebra of rank 2n 

{T IM ,T V } = 26^ (113) 
and they can then be used to construct the generators of 50(2n) as follows: 

s^ = i [r^r,]. (ii4) 

The dimensionality of the spinor representation of 50(2n) is 2™. In terms of the 
SU(n) basis, the spinor representation of 50(2n) can then be constructed by, 

|0>~1 (115) 

|0 > ~ n (116) 

eJ ej io > - n -^p± (H7) 
el e] ^ io > ~ "("-y"- 2 ) (us) 

(H9) 

^ ... & |0 > ~ n (120) 
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where |0 > is the SU (n) invariant vacuum state. The spinor representation can then 
be split into two 2™~ ^dimensional representations by a chiral projection operator. 
Let us define 

r = »Ti r 2 ... r 2n (121) 

and the number operator 

Ni = &i. (122) 

r then can be written as 

r = 

n 

= n (! - 2n *) 

i-l 

= (-1)™. (123) 

To arrive at the last step, we have used the property of the number operator nf — 
Hi to get 1 — 2n,i — (— \) ni and n — ^ rij. One can then check that 

P^,r o ]=0. (124) 

The chirality projection operator is therefore defined by 

^(i±r ). (125) 

Consider the case n = 5 and define a column vector \ip >: 

|V> > = |0 > + <t|0 > iPj + l&t\0 > iP jk + ^e ij ' Wm eki + 4|0 > 

+ le*'™&t&&|0 > ^ + ^t^t| Q > ^ (126) 



where ^ is not the complex conjugate of ip but an independent vector. This can be 
generalized to any n if we write 

^=(^^i^^?i?o) T ' ( 127 ) 



The spinor representation is then split under the chirality projection operator as 



* = ( ! + ) (128) 



where 



and 



^± = l(i±r )v (129) 



= *i , 1>- = • (130) 
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In the case of n = 5, ip i and ipij are 5 and 10-dimcnsional representations of 5(7(5) 
and tpo is the singlet. All the SM fermions are assigned to The electric charge 
formula for 50(10) is given by 

Q= Is 7 l(£ 12 + £ 34 + £ 56 ). (131) 
2 o 

The dimensionality of the adjoint representation of SO(2n) is ( 2 ")(^"~ 1 ) _ p or 
50(10), it is 45-dimensional. Under the decomposition with respect to 5(7(3) x 
SU(2) L x SU{2) R these 45 gauge bosons are: 

45 = (8, 1,1) + (1,3,1) + (1,1,1) 

+ (1, 1, 3) + (3, 2, 2) + (3, 2, 2) + (3, 1, 1) + (3, 1, 1). (132) 

In this basis, the 12 Standard Model gauge fields are in the (8, 1, 1), (1, 3, 1) and 
(1, 1, 1) multiplets. The rest are 33 new gauge bosons which could potentially me- 
diate proton decay. 



4.2. Symmetry Breaking 

Because 5O(10) is a rank-5 group while SM is a rank-4 group, there exist several 
intermediate symmetries through which 50(10) can descend to 5(7(3) x SU(2)l x 
U(1) Y - There are four maximal subgroups of 5O(10): 5(7(4) x SU(2) L x SU(2) R , 
SU{5) x U(l), 50(9), and 50(7) x 5(7(2). Only through the breaking chains of 
5(7(4) x SU(2) L x SU(2) R and 5(7(5) x (7(1) can one obtain the correct quantum 
numbers for the SM particle content. Due to the presence of these intermediate 
scales, the predictions for proton lifetime and sin 2 9 W are much less certain, com- 
pared to the case of 5(7(5). Details of breaking chains giving the SM are as follows: 
(i) The left-right symmetry breaking chain is 

50(10) <^> 5t/(4) x SU(2) L x SU(2) R 

<4 ^ 5(7(3) x SU(2) L x SU(2) R x U(1) B -l 
<i26ei26> SJJ ^ x SU ( 2 ) L X U(1) Y 

<^> 5(7(3) x U(1) EM - (133) 

In this case, the hypercharge Y is given by Y — 2T 3R + (B — L). The first step 
of breaking to the left-right symmetry group is achieved by a symmetric two-index 
tensor, < 54 >. 5(7(4) in the left-right symmetry group is then broken to 5C/(3) x 
U(1)b-l by the adjoint < 45 >. The subsequent breaking to the SM gauge group 
is achieved by the anti-symmetric 5-index tensor, < 126 > and < 126 >; the 
electroweak symmetry breaking is then achieved by < 10 >. In realistic models, 
the two Higgs doublets in MSSM are linear combinations of the 5(7(2) doublet 
components from different 50(10) representations of Higgses, for example, 10 and 
126. Thus all fields in the linear combinations contribute to electroweak symmetry 
breaking. 
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(ii) For the SU (5) breaking chain, 

50(10) -» SU{b) x 1/(1),. 

-> 517(3) x 5*7(2) x C7(l) z x J7(l) x 
-> 517(3) x 5*7(2) x J7(l) y 

-> 517(3) x U{\) EM (134) 

the hypercharge Y is given by ^Y = az + /3x where z and x are the charges under 
U(l) z and f7(l) x respectively. There are two possible ways to embed the SM under 
this route: (a,/3) = (1/6,0) or (-1/15,-2/5). In the case of (a,/3) = (1/6,0), we 
obtain Y — \{2z). This corresponds to the 5t7(5) breaking chain 

5O(10) < 16 ® T 6> SU(5) 

<4 ^> 517(3) x SU(2) x U(1) Y 

<^> SU{3)xU{l) EM . (135) 

In this case, the spinors < 16 ® 16 > break the symmetry down to 5*7(5); < 45 > 
then breaks SU(5) down to the SM. In the case of (a, 5) = (—1/15,-2/5), we 
have Y = =±(z + &x). This corresponds to flipped 5f7(5) (that is, 517(5) x U(l)) 
breaking chain 

5O(10) 517(5) x U(l) x 

<i 6 el6> su{i) x SU ( 2 ) x jj^ y 

<^> 5*7(3) x U(1) EM - (136) 

For this breaking to occur, again we need < 16 16 > and < 45 >. Each breaking 
chain has its characteristic set of Higgs fields and symmetry breaking superpoten- 
tial. The electroweak symmetry breaking is then achieved by < 10 >. In realistic 
models, the two Higgs doublets in MSSM are linear combinations of the SU(2) 
doublet components from different 5O(10) representations of Higgses, for example, 
10 and 16. As we will see below, different symmetry breaking chains give rise to 
different mass relations between various sectors. 



4.3. Renormalization Group Equation and the Georgi-Jarlskog 
( GJ) relations 

Before describing the Yukawa sector of the 5O(10) models, we discuss how to relate 
the weak scale observables to the GUT scale p arameters. We use the expressions 
derived from 1-loop RGEs given byESEZESESI 

m u = Y^RuijuB^Vu, m c = Y^R u i] c Bfv u , m t = Y°R u B^v u 
m d = Y*j>R d r) d v d , m s = Y°R d rj s v d , m b = Y b R d r] b B t v d 

m e =Y°R e v d , m^ = Y°R e v d , m T = Y?R e v d (137) 
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and 

V. ij = ud,us,cd,cs,tb 

Vij ~ \ V? 3 Bt \ ij = ub, cb, td, ts ( M) 

where Vij are CKM matrix elements; quantities with superscript are evaluated at 
GUT scale, and all the to/ and Vij are the experimental values at Mz lUZESEll yo 
are Yukawa couplings, and v u and v d are VEV of the Higgs fields H u and H d . The 
SM Higgs VEV is v — \/ v u + v d = 246/V2 GeV. The running factor 77/ includes 
QCD + QED contributions: For / = b,c, rjj is for the range to/ to m t , and for 
/ = u,d, s, rjf is for the range lGeV to m t ; 

Vu = Vd = Vs = 2.38t°;^ 

Vb = 1.53lg;g. (139) 

Ru,d,e are contributions of the gauge-coupling constants running from weak scale 
M z to the SUSY breaking scale, taken to be mt, with the SM spectrum, and from 
to 4 to the GUT scale with MSSM spectrum; 

^ = 3.5318:8?, R d = 3.4318:85, ^ = 1.50. (140) 

B t is the running induced by large top-quark Yukawa coupling defined by 

1 



B t = exp 



16tt 2 



In Mgut 

Y^Mln/i) 

In Ms crsv 



(141) 



which varies from 0.7 to 0.9 corresponding to the perturbative limit Y t ° w 3 and 
the lower limit Y t ° rj 0.5 imposed by the top-pole mass. 

Naively, one would expect that masses of the down type quarks are identical to 
masses of the charged leptons, because of the quark-lepton unification. This turns 
out to be not true. Taking the experimentally measured values for masses of the 
down-type quarks and those of the charged leptons, evolving these values using 
the RGE's from the weak sc ale to the GUT scale, however, one finds the following 
approximate relationa22E2l 

m d ~ 3to 6 , m s ~ ^to m , m b ~ m T . (142) 

These are known as the Georgi-Jarlskog relations. As it will become apparent later, 
one way to satisfy these relations is by having a relative factor of —3 in the (22) entry 
of the charged lepton mass matrix with respect to that of the down-type quarks, and 
all other elements are identical in these two mass matrices. An example suggested 
by Georgia and Jarlskog is: 

/0 A 0\ / E 0\ /0 E 0\ 

M u = I A B , M d = I E F , M e = I E -3F . (143) 
\0 BCj \0 GJ \0 Gj 

This factor arises naturally as a Clebsch-Gordon (CG) coefficient in some models 
of 5*0(10). 
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4.4. Yukawa sector in 50(10) 

In 50(10), at the renormalizable level, only three types of Higgs fields can couple 
to fermions, 



16® 16 = 10 s © 120 A © 126s 



(144) 



namely, 10, 120^, and 126s, where the subscripts 5 and A refer to the symmetry 
property under interchanging two family indices in the Yukawa couplings y a b. That 
is, 



126 



-i; 10 _ -nio v 120 = —V 120 V 126 = V 

•s ab •Sba i ^ ab ^ba i ^ ab ^ba 

The gauge invariant Yukawa couplings are then given by 

y a 1 °(16) a (16) fc (10) = y^ T a BC^Y a ^ b H a 
y a 120 (16) a (16) b (120) = y^ T a BC- x Y a Y (i Y^ h k a(il 

3^ 26 (i6) a (i6) 6 (T26) = y a 1 6 26 *rso- 1 r a r /3 r 7 r 5 r l * 6 A a/ 3 7 ^ 



(145) 

(146) 
(147) 
(148) 



where 'J denotes the matter fields; H ai A Q ^ 7 and A 0| g 7 ^ denote the 10-, 120- and 
126-dim Higgs fields, respectively. C is the usual Dirac charge conjugate operator 
and B is the charge conjugate operator for 5O(10) defined by 



and we can choose 



B : 



n 



(149) 



(150) 



=odd 



Under this charge conjugation 



B 



{ TP0\ 




( A ) 














Ipi 

i>ij 










\ ^ ) 



(151) 



Note that 5O(10) can break down to SM through many different breaking 
chains. Different breaking chains give rise to different mass relations among the 
up-quark, down-quark, charged lepton and neutrino sectors. In what follows, we 
discuss the two symmetry breaking separately. 

(i) SU(5) breaking chain: 

Under the SU(5) decomposition, we have 

16=1 + 5 + 10 (152) 
10 = 5 + 5 (153) 
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120 = 5 + 5 + 10+ 10 + 45 + 45 (154) 
126 = 1 + 5 + 10+ 15 + 45 + 50 (155) 

where the SU(5) component 5 and 45 contain the Y = +1 SU{2)l Higgs doublet, 
and the SU(5) component 5 and 45 contain the Y = — 1 SU(2)l Higgs doublet. 
The SU(5) singlet contained in 126 gives masses to the RH neutrinos through 
the coupling (16j)(16j)(126ij). (As we will discuss later, some models utilize non- 
renormalizable operators (16 i )(16 J )(16#)(16H) to generate RH neutrino masses. 
This can be achieved because 16h also contains a SU(5) singlet component.) The 15 
of SU(5) contained in 126 has a (1, 3, 2) component under SU(3) x SU(2) L x U(l)y 
which couples to two lepton doublets as (1, 2, — 1)(1, 2, — 1)(1, 3, 2)h under Gsm and 
gives the LH neutrino Majorana masses in the Type II see-saw mechanism. 

As the neutral components in these SU(2)l doublets acquires VEV's of the 
electroweak scale, the following mass matrices are obtained 

M u = yll (5(10)) + y«o (45(120)) + y™ (5(126)) = Y u v u (156) 
M d = yH (5(10)) + ((5(120)) + (45(120))) + (45(126)) = Y d v d (157) 
Me = yll (5(10)) + ^ 6 20 ((5(120)) - 3 (45(120))) - 3^(45(126)) = Y e v d (158) 
M VLR = yll (5(10)) + y™ (5(120)) - 3^If (5(126)) = Y VLR v u (159) 

where we denote the m-dim SU (5) component of the n-dim representation of 
5O(10) by m(n); v u and v d are the vacuum expectation values of the two Higgs 
doublets in MSSM. A Clebsch-Gordon coefficient (—3) is generated in the lepton 
sectors when the SU(5) component 45 are involved in the Yukawa couplings. This 
factor of (—3) is very crucial for obtaining the Georgi-Jarlskog relations as we have 
seen in the previous section. The neutrino Majorana mass matrices are given by 

M ViRR = y^ (1(126)) (160) 
M VtLL =y\f (15(126)). (161) 

To see how these CG coefficients (—3) come about, let us decompose the fol- 
lowing representations under SU(5). The 120 and 126-dimensional representation 
of 50(10) both contain a component which transforms as 45 under SU(5). The 
Yukawa interactions then can be written as 

(16)(16)(120), (16)(16)(I26) D (10) (5) (45) = V Q ^ 7 <Kp>- (162) 
We then write out all the terms in the summation 

V> a/3 V7 < Kfi >^ (V' 45 ^ < ^45 > +r b i>a < H a ab >) = (-3e+e- + dd c )v . (163) 

Here the index a = 1,2,3. Note that H 1 ^ is anti-symmetric under inter-changing 
a <-» (i and it is traceless 

H2 = -Hj a , ( 164 ) 
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This implies that 

3# Q a 5 + #45 = 0, < Hi, >= -3 < H a a5 -3v' . (165) 
Essentially, the CG factor of (—3) is related to the fact that there are three colors. 

We note that if 10 and 126 are the only fields utilized in the Yukawa sector, we 
have the up-quark mass matrix related to the Dirac neutrino mass matrix, and the 
down-quark mass matrix related to the charged lepton mass matrix. When 120 is 
introduced, the relation between the up-quark sector and the Dirac neutrino sector 
is lost because these two sectors receive contributions from different components of 
120. 

(ii) SU(4) x SU{2) L x SU(2) R breaking chain: 

Under SU(4) x SU(2) L x SU(2) R , the relevant 5O(10) representations have the 
following decomposition 

16 = (4, 2,1) + (4, 1,2) (166) 

10= (6, 1,1) + (1,2, 2) (167) 

120 = (15, 2, 2) + (6, 3, 1) + (6, 1, 3) + (1, 2, 2) + (10, 1, 1) + (TO, 1, 1) (168) 

126 = (10, 1, 3) + (10, 3, 1) + (15, 2, 2) + (6, 1, 1) (169) 

where the components (15,2,2) and (1,2,2) both contain a pair of the Y = ±1 
SU (2) l Higgs doublets, whose neutral components give masses to the fermions. The 
component (10, 1, 3) contained in 126 gives masses to the RH neutrinos through the 
coupling (16j)(16j)(126jf). The LH neutrino Majorana masses are generated due 
to the (10,3,1) component of 126. As these SU(2)l doublets acquire VEV's, the 
following mass matrices are generated, 

M u = y% (10+) +^ b 20 «120 + ) + I (120'+)) + \y™ (l26 + ) = Y u v u (170) 

M d = yll (io-) + y%°{- (i2o-) + 1 (120'-)) - \y^{m~) = v d v d (171) 
Me = yll (10-) + y a l°{- (120-) - (120'-)) + y}f (m~) = Y e v d (172) 

M VLR = yll (10+) +^f((l20+) - (l20'+)) +y™ (T26 + ) = Y VLR v u . (173) 

Note that a Clebsch-Gordon coefficient (—3) is generated in the lepton sectors when 
the SU(4) x SU(2)l x SU(2)r components (15,2,2) are involved in the Yukawa 
couplings. The neutrino Majorana mass matrices are given by 

M v ,RR=y^ (T26~'°) (174) 

M v , LL =y~T (l26 ,+ ) (175) 
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Table 3. Standard Model singlet components contained in various 50(10) Higgs rep- 
resentations. Here the subscripts refer to the symmetry groups on the right-hand 
side of Eg. 11331 : and superscripts +/0/— referring to the sign of the hypercharge Y. 



(10+) 


(1,0)31 


c 


(1,2, 1)321 


C (1,2,2,0)3221 


C (1,2,2)422 


C 10 


(io-) 


(1,0)31 


c 


(1,2,-1)321 


C (1,2,2,0)3221 


C (1,2,2)422 


C 10 


(120+) 


(1,0)31 


c 


(1,2, 1)321 


C (1,2,2,0)3221 


C (1,2,2)422 


C 120 


(120-) 


(1,0)31 


c 


(1,2,-1)321 


C (1,2,2,0)3221 


C (1,2,2)422 


C 120 


/l20'+\ 


(1,0)31 


c 


(1,2, 1)321 


C (1,2,2,0)3221 


C (15,2,2)422 


C 120 


(l20'-\ 


(1,0)31 


c 


(1,2,-1)321 


C (1,2,2,0)3221 


C (15,2,2)422 


C 120 


/l26 + \ 


(1,0)31 


c 


(1,2, 1)321 


C (1,2,2,0)3221 


C (15,2,2)422 


C 126 


(T26-) 


(1,0)31 


c 


(1,2,-1)321 


C (1,2,2,0)3221 


C (15,2,2)422 


C 126 


(T26°) 


(1,0)31 


c 


(1,1,0)321 


C (1,1,3,-2)3221 


C (10,1,3)422 


C 126 


(T26 ,+ ) 


(1,0)31 


c 


(1,3,2) 321 


C (l,3,l,2)3 22 i 


C (TO, 3, 1)422 


C 126 



where various VEVs are those of the neutral components of 5*0(10) representations 
as indicated in Table |3J 

Non-renormalizable Operators: 

If 5*0(10) breaks down through the 5/7(5) breaking chain, the Higgs fields 
needed are in the (16 © 16), 45 and 54-dim representaions. A lop-sided texture 
can arise, when non-renormalizable operators involing the 16's are utilized to gen- 
erate fermion masses. The 16-dim Higgs fields are needed in this case to break 
5O(10) down to 5/7(5). A lop-sided texture is generated by the operatoJ^ 

A(16 i 16 Hl )(16 i 16 H2 ). (176) 

If 16h± acquires a VEV along the 5/7(5) singlet direction which breaks 5O(10) 
down to 5/7(5), and 16^ 2 acquires a VEV along the 5 direction of 5/7(5) which 
breaks the electroweak symmetry 5/7(2) x /7(1), we obtain 

(5i)(10j) <l Hl ><5 Hl >■ (177) 

Inside the first parenthesis in Ea. (jl76JI . the two 16's contract to from a 5 of 5/7(5), 
while inside the second parenthesis, the two 16's contract to form a 5 of 5/7(5). 
This contraction arises by integrating out a pair of 5 and 5 of 5/7(5) from the 10's 
of 5O(10), as shown in Fig'QJ Because the 5 contains the 5/7(2) lepton doublet and 
the singlet down-type quarks, the resulting mass terms 

Kdl,id L ,i + e L ,iel d )v d (178) 

are related by Md = Mj and the lop-sided mass texture arises, 
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if is chosen to be (2, 3). The (33) entry which is expected to be of order 0(1) 
is generated by tree level diagram involving a 10. The (23) and (32) entries of 
Md and M e also receive contributions from other non-renormalizable operators, for 
example, 16^16^45^10//. If A is of order 0(1) and other contributions to the (23) 
and (32) entries are much smaller than one, a large mixing angle is then generated 
in the right-handed down quark sector, and the left-handed charged lepton sector, 
while the corresponding mixing angle in the left-handed down quark sector is small. 
This thus provides a way to explain the large mixing angle in atmospheric neutrinos 
while the quark mixing V c b is small. 

As we have seen above, when additional matter fields are introduced into the 
model and non-renormalizable operators that generate fermion masses are taken 
into account, as in the case of the Froggatt-Nielsen mechanism, other Higgs rep- 
resentations can play a role in the mass generation. An interesting case is the 
45-dimensional Higgs representation which has Yukawa couplings to a 16- and a 
10-dim ma tter fields. These non-renormalizable operators can be expressed generi- 
cally bjEO 



O tj = 16, • O n • 10 • m • 16, 
where the operator 0„ is given by 

■45„ 



= Mg45 p+1 . 



n, p, q — integer. 



(180) 



(181) 



Here Mq and Mpi refer to the GUT scale and the Planck scale respectively. The 
45-dimensional representation Higgs can acquire VEV along the following four di- 
rections: X, Y, B — L, T^ji, where X and Y are defined as 



X = -(Eia- 
Y = y<g>C 

where 



^34 



c 



^56 T 

-% 

1 



^78 



^9 10 



(182) 



y = diag(l/3, 1/3, 1/3, -1/2, -1/2). (183) 



A systematical way to search for these effective operators is discussed in Ref. 95. 
The coefficient —3 is obtained in this case when the 45-dimensional Higgs acquires 



l(16 Hl ) 1(1h) 5(16hJ 



5(160 5(10) 5(10) 10(16j) 

Fig. 4. Froggatte-Nielsen diagram which generates lop-sided texture. Here we denote the field 
by m(n) where m is the SU (5) representation which is a component of the n-dim representation 
of SO(10). 
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a VEV along the B — L direction. We summarize in Table. 0] all possible CG 
coefficients for the four possible directions of < 45 >. 



Table 4. CG coefficients for the four possible directions of 
< 45 >. 





a 


u 


d 


d 


e 


e 


V 


;/ 


B-L 


1 


-1 


1 


-1 


-3 


3 


-3 


3 







-1/2 





1/2 





1/2 





-1/2 


X 


1 


1 


1 


-3 


-3 


1 


-3 


5 


Y 


1/3 


-4/3 


1/3 


2/3 


-1 


2 


-1 






4.5. Automatic R-parity conservation and 126 v.s. (16 X 16) for 
neutrino masses 

One of the salient features of 5*0(10) is that in some classes of models, R-parity 
is conserved automatically. Group theoretically, the congruence number is defined 
as follows: For an irreducible representation of SO(2n) whose Dynkin index reads 
(ai a 2 ...a n ), this representation has congruence number 

(ci, c 2 ) = (a„_i + a„, 2a x + 2a 3 + ... + 2a„_ 2 + (n- 2)a„_ 1 + na n ) 

mod (2,4) for n= odd. (184) 

There are four possible classes of (01,02): (0, 0), (0, 2), (1, 1), (1, 3). The first two 
classes are tensor-like while the later two classes are spinor-like which is trouble- 
some. To see how (ci, c 2 ) relates to R-parity, alternatively we can define, in S'O(IO), 

c = 3(5 - L) mod 4 (185) 

and c = Max (01,02). It has been shown that if all the Higgs representations that 
break 50(10) down to the SM are chosen to have congruence number c = or 2, 
then R-parity is preserved at all energiesP^^I Representations having c — are: 
45, 54, 210,...; those having c = 2 are: 10, 126, 126,.... Note that the spinor 
representations 16 and 16 have c = 1 and 3 respectively. 

Some models avoid the use of 126-dim Higgses by introducing non- 
renormalizable operators of the form jg-?/'a'0&(16.ff )(16h ) instead of a renormalizable 
ipa^b^^H- Models utilizing the spinor representation 16 to construct the neutrino 
mass operators generally have R-parity broken at some high energy scale. Such mod- 
els also appear to be less constrained due to the inclusion of non-renormalizable 
operators. Also, a discrete symmetry, the R-parity symmetry, must be imposed 
by hand to avoid dangerous dim-4 baryon number violating operators in the ef- 
fective potential at low energies which otherwise could lead to fast proton decay 
rate. Utilizing 126-dim representation of Higgses has the advantage that R-parity 
symm etry is automatic . The 1 2 6 repr esent at ion has been used in model building 
before! 98 l 9iJ l 1UU l 1U1 l 49 l iU2 l 48 l 1U3 l 1U4 llt is to be noted that the contribution of the 
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126-dimensional representation to the /3-function makes the model nonperturbative 
(with the onset of the Landau pole) above the unification scale Mqut- One could 
view these models as an effective theories valid below this scale where coupling 
constants are perturbative. One may also argue against the use of 126 with the 
fact that it is not possible to obtain such a large representation from string theory. 
It has been shown that in heter otic s tring theory it is not possible to get 126 of 
50(10) up to Kac-Moody level-5^^ Nevertheless, no such constraints have been 
found in other types of string theoriesP23 



4.6. Some Related Issues 
4.6.1. Proton Decay 

One of the signatures of any GUT model is baryon number violating processes. 
These processes include proton decay (A(B—L) = 0), N—N oscillation (A(B—L) ^ 
0), etc. The theoretical prediction for oscillation time in N — N oscillation can natu- 
rally satisfy the experimental lower limit T N _j7 > 0.86 x 10 8 secpHlif a high B — L 
breaking scale is assumed. On the other hand, the non-observation of proton decay 
has put many GUT models under siege. There are three kinds of operators leading 
to proton decays in SUSY GUT's: 

(i) Dimension-6 operators: 

As we have mentioned previously, the extra gauge bosons in GUT models, the 
X and Y gauge bosons, can lead to proton decay. The terms in the Lagrangian 
containing the X Y gauge bosons are 

£ x ,y = -^X^i(eijkUk L 7iJ,UjL +d i j t _ l e + ) 

+ l -^Y t , 4 (e ljk u c kL -/ tl d jL - Uii7 M eJ + d lR ^^v c R ) + h.c. (186) 

and the following vertices are allowed 

~X^uu,~de + , Y->ud,ue + . (187) 

These can thus lead to proton decay via the dim-6 operators. Note that these type 
of operators exist in both non-SUSY and SUSY GUT's. The domin ant d ecay mode 
is p — ► e + 7r°, and the decay amplitude associated with this mode i d^^ l 

M p ^ e +^o ~ 4kolgvt/M% ut (188) 

leading to a life-time off^EI 

^ 1 a g v -|n29 + 0.7; M GUT ^ 4 hRQ s 

P MU^ P " ( 2.1xlO"GeV ) (189) 

where m p is the mass of the proton. For Mqut — 2 x 10 16 GeU, we get r p ~ 
4.5 x io 37 - ' 7 years which is far above the current capability of SuperKamiokande 
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experiments whose limit is ~ 10 years. 

(ii) Dimension-5 operators 

In SUSY GUT's, a new channel for proton decay is possible via the dim-5 operator 
through the exchange of the color triplet Higgsinos where QQH and QLH via HH 
mixing generate an effective operator 

QQQL/Mh. (190) 

In order to suppress this operator, the masses of the color triplet Higgsinos must 
have superheavy masses. From the point of view of unification, we would like to 
have the spectrum of MSSM below Mqut- This requires that all the Higgs fields, 
including the color triplet Higgsinos, to be very heavy, with masses of the order 
of Mqut, except the pair of SU(2)l doublet Higgses remaining light which are 
then identified as the pair of Higgs doublets in MSSM. How to achieve such a mass 
splitting is referred to as the doublet-triplet splitting (DTS) and doublet- doublet 
splitting (DDS) problems. The dominant de cay mode of these operator is r(p — > 
K + v) ~ m~ ~ Mq Ut . Its decay amplitude iiEPJ 

\A h u hd m g augino9GUT 

M *"« + »- M B IG**M% ' (191) 

Dimopolous and Wilczek proposed a mechanis m^ ^ 1 to achieve such mass splittings 
using < 45 h > along the (B — L) direction 

< 45 H >— ir 2 ®diag(a, a,a,0, 0). (192) 

Chacko and Mohapatr JH^ found that with a complimentary VEV pattern to that 
of the Dimopoulos- Wilczek type, that is, 

<45 H >=iT 2 <8> diag(Q,0,0,b,b) (193) 

the same goals can be achieved. A solution utilizing 126 h instead of 45 h to achieve 
the DTS and DDS was proposed by Lee and MohapatraPSI Detailed calculations 
have shown, nevertheless, that even with these mechanisms in place, in order for 
the prediction of t p to be consistent with the experimental limit, the effective Mh 
must be larger than Mqut by at least a factor of 10. This in turn requires some 
couplings to be much larger than 1 which is somewhat unnatural. 

(iii) Dimcnsion-4 operators: 

As we have seen, the dim-4 operators are forbidden if there is R-parity in the 
model. In SUSY SU(5) one has to impose R-parity by hand, while in some S'O(IO) 
models, R-parity is automatically conserved if certain type of Higgs fields are chosen 
to construct the model, as we have discussed in the previous section. 
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4.6.2. Baryogenesis 

The three Sakharov conditionsj ^^ (i) baryon number (B) violating processes, (ii) 
C and CP violation, and (Hi) the deviation from thermal equilibrium, for the gener- 
ation of the cosmological matter anti-matter asymmetry can be naturally satisfied 
in the SO(10) model. Due to the presence of sphaleron effects, the only chance for 
GUT baryogenesis to work is to produce an asymmetry in B — L at a high scale. 
To see this, let us first write the baryon number B as 

B = ±{B + L) + ~(B-L). (194) 

The electroweak sphaleron transitions rapidly erase the asymmetry B + L as soon as 
the temperature drops down to about 10 12 GeV. Therefore, to have a non- vanishing 
baryonic asymmetry requires a non- vanishing B—L asymmetry. The crucial point to 
note is that even though the B+L asymmetry is erased by the sphaleron transitions, 
the orthogonal combination B — L is left untouched, and it opens up the possibility 
that the baryonic asymmetry can be generated through leptogenesis F^^^^^ 
The basic idea is that, since B + L must vanish at all times due to the sphaleron 
transitions, the asymmetry in the lepton number will consequently be converted 
into the asymmetry in the baryon number (with an opposite sign) . The primordial 
leptonic asymmetry is generated by the out-of-equilibrium decay of the heavy right- 
handed Majorana neutrinos and their scalar partners in the supersymmetric case. 
The relevant superpotential is 

1 

2 

The heavy right-handed neutrinos and their scalar partners can decay through the 
following four decay modes: 



= (Y e )ijEfLjHi + (Y UL r ) ij N?Lj H 2 + - (M RR )ij N?Nj. (195) 



N x — ► I + ti c (196) 

Ni^l + H 2 (197) 

N£ — >l + H 2 (198) 

Nf — > I + h c . (199) 

The interference between the tree-level and one-loop diagrams thus gives rise to the 
CP asymmetry. 

In the basis where both charged lepton Yukawa couplings and the right-handed 
neutrino mass matrix are diagonal, the amount of CP asymmetry due to the inter- 
ference between the tree level and one-loop diagrams for each decay mode is given 

bJ imi 

1 1 Ju[2 

»7T (h v hl)u - 2 , M l 



where 



/(x) = VS[M— ) + r^r]. ( 201 ) 
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and h v is the Dirac neutrino Yukawa matrix in the new basis. The right-handed 
Majorana neutrino mass matrix is diagonalized by 

PmOrM rr O t r P m = diag{M u M 2 , M 3 ) (202) 

where Mi's are real and non-negative, and Pm is the diagonal Majorana phase 
matrix. In terms of the diagonalization matrices and the original Dirac neutrino 
Yukawa coupling, (h v h\) can then be rewritten as 

Khl = P MRR R Y VLR Ut L UeA R R P^ RR = Pm rr O r Y Vlr Y} lr O r P Mrr . (203) 

We see that the phases in the right-handed neutrino mass matrix M RR and the 
Dirac neutrino mass matrix Y VLR are responsible for the CP asymmetry needed for 
the leptogcnesis. For a hierarchical heavy right-handed neutrino mass spectrum, 
M3 3> Mi ^ Mi, the argument of the function f(x), x = is much greater than 
1. We can then approximate f(x) as 

/W = ^[(i-i + -) + ^T]-^ = 3 t- (204) 

The asymmetry t\ can thus be rewritten as 

ei * -^TTIJV- E Im {{h v hlf u }{^) (205) 
»7T (h„hl)u ,f^ 3 M i 

where the factor 4 accounts for the fact that there are four decay modes. Using the 
fact that the mixing in Y VLR is small and that (Y„ LR ) 33 dominates other elements, 
we can further approximate 

(h v hl) u c (PnPu 1 )^ \ 2 (0 R )i 3 (0* R ) i3 (206) 
(Khl) n ~ |y„ T | 2 |(O fl )i 3 | 2 (207) 

and 

— — — r- ^\y VT \ lm{{F n F u ) — (0 R ) i3 \. (208) 

(h v hl)u (Or) 13 

To have a large amount of CP asymmetry, ei, thus requires that the hierarchy 
among the three right-handed neutrino masses cannot be too large (that is, 

cannot be too small), and that the imaginary part of {(fii-P^ 1 ) 2 {o R y 3 (O r ) 2 3 } 
together with the neutrino Dirac Yukawa couplings cannot be too small. 
The amount of the lepton asymmetry generated is given by 
y L ^!^n lL = K ei 
s 5* 

Here g» is the number of relativistic degrees of freedom. For MSSM, it is = 
228.75. The out-of-cquilibrium decay of the heavy Majorana neutrinos requires the 
decay width of the lightest neutrino, Ti, smaller than the Hubble constant at the 
temperature of the decay. That is, 

ri M p i (h v ht)u 



H\ T=Ml (1.7)(32tt)V^ M 1 



< 1 (210) 
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where Mpi is the Planck scale taken to be 1.2 x 10 19 GeV. In general, one can still 
have a sizable CP asymmetry remains even for 1 < r < 10. The wash-out effects 
due to inverse decays and lepton number violating scattering processes together 
with the time evolution of the system is then accounted for by the factor n. It is 
obtaine d by solving the Boltzmann equation for the system. An approximation is 
given b y I 

10 6 <r: k = (O.lrf/Se-CfXo.ir) 1 ^ (< HT 7 ) (211) 
10<r<10 6 : K =-^ JSX (10- 2 -10- 7 ) (212) 

0<r<10: K =57=ITI (lO" 1 - lO" 2 ) (213) 

where inside the parentheses we give the order of magnitude of k for each cor- 
responding r. We note that in order to have a small dilution factor, the lightest 
right-handed neutrino cannot be too light. The electroweak sphaleron effect will 
convert the lepton asymmetry Yj, into baryon asymmetry Yg , and they are related 

by 

Yb = = cYb—l = ^~Y L (214) 

s c — 1 

with 

8N F + AN H 



(215) 
22iV F + 13N H 

where Np is the number of families and Njj is the number of SU (2) Higgs doublets. 
For the MSSM spectrum, (Np, Nh) = (3, 2), we have the conversion factor (^zt) — 
-0.53. 

Models with only two RH neutrinos have been constructed by Frampton et al&^ 
and by Raby^l (see Sec. 3 and 6), which give rise to bi-large mixing pattern and a 
correlation between the sign of the baryon number asymmetry and the sign of the 
CP violation in neutrino oscillation. General analyses on the consistency between 
constraints from simplest 50(10) models and leptogenesis can be found in Ref. 117, 
118. 



4.7. SUSY Breaking 

SUSY breaking can be incorporated into models by including explicitly the soft 
SUSY breaking terms. Since the RGE's for the Yukawa coupling constants and 
gauge coupling constants do not have any dependence on the soft breaking param- 
eters up to two-loop level, the presence of these soft SUSY breaking terms does not 
affect the predictions for fermion masses and mixing angles. On the other hand, 
since the evolutions of soft SUSY breaking parameters does depend on the Yukawa 
coupling constants and gauge coupling constants, whether the EW symmetry is 
broken (that is, the mass-square of the light Higgs doublet is driven to be negative) 
may depend on the Yukawa sector, which can be used as a test of the validity of 
the model. 
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H* TT 

u M u 

X X 

I I 

I I 

I I 

1 i 



^g^gjg^pl t„=<g> 



Fig. 5. Leading order contributions to the threshold corrections to the mass of the b-quark. 

Even though SUSY breaking does not affect the running of the gauge coupling 
constants and that of the Yukawa coupling constants, it could have a large contri- 
bution to the threshold corrections, which is the subject of the next section. 



4.8. Threshold Corrections 

When the RGE analysis is performed, we usually consider (3 function coefficient 
as a constant for each coupling constant between the two relevant scales. This is 
done under the assumption that all the heavy modes decouple at the same scale, 
the symmetry breaking scale. Threshold corrections are the corrections due to the 
differences between the symmetry breaking scale and the masses of the heavy par- 
ticles decoupled from the spectrum after symmetry breaking takes place. If all the 
heavy particles acquire masses exactly the same as the symmetry breaking scale, 
there are no threshold corrections. In practice, this is not the case. There are two 
possible sources of threshold corrections: 

(i) GUT scale threshold corrections: Due to the presence of many large Higgs rep- 
resentations in 50(10), the GUT scale threshold corrections could be large. 

(ii) SUSY threshold corrections: Large threshold corrections to m&, Srrib/mb ~ 
— (0.15 ~ 0.2), are needed in most 50(10) models in order to have a prediction for 
TOh consistent with experiment. The dominant contributions are from the diagrams 
shown in Fig. They give a correction 

8m b I m b ~ (tan/3/50)/ (216) 

and / is given byl 119 l 12U l 12 H 

' ~ 50 UG r ° ^ 9if(^) - 2A?/(-£-)] (217) 



167r 2 a2 miff 3c*g' 



m eff ^ 1 m eff 



where f(x) = (1 — x+xlnx)/(l — x) 2 ; trig and \x are the gluino mass and the [i term 
evaluated at the weak scale.; m^ff = \{ m i + m |,) is the average of the squared 
masses of the 5f7(2)-singlet bottom squark and the 5t/(2)-doublet third genera- 
tion squarks. A large soft SUSY breaking parameter space can give rise to such a 
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correction. With the typical values aa — 0.75 and (a 3 ,a 2 ,y 33 ) — (0.124,0.034,1) 
at M w 

eaki $b can then be approximated as, with the assumption 2777-1/2 — ^o^^W' 

I = (0.315) x t (0.69/(t 2 ) - f(x 2 )) (218) 

where t = (///mo) and x = m 1 / 2 /mo with m 1 / 2 and mo being the gluino mass and 
the common scalar mass respectively in the constra ined MSSM (CMSSM). With 
(t,x) — (5.2,2) which are typical values in cmssmEU a correction (Snib/mi,) = 
—0.15 for tan/5 = 10 is obtained. 

5. SO(10) Models with Texture Assumptions 

In what follows we discuss a few selected £0(10) models combined with some 
texture ansatz for the mass matrices. 



5.1. Buchmuller and Wyler 

Buchmuller and Wylei^^ assume symmetric mass textures for the up- and down- 
type quarks 

/ eV* \ 
Mu,d ~ e 3 e^ pe 2 rje 2 . (219) 
V r,e 2 ) 

Here p = \p\e la and r/ = \r]\e 1 ^ are complex parameters of 0(1). Using the 50(10) 
relations, they have 

m^rac = m ^ ^ = md (220) 

assuming the incorrect mass relations in the lighter two generations are lifted when 
higher dimensional Higgs representations are introduced. The right-handed neutrino 
Majorana mass matrix is generated by 126//, and is assumed to have the following 
form 

/ Mi 2 \ 
M UtRR = Mi 2 M 22 M 23 . (221) 
\ M 23 M 33 ) 

With the relations 

M12 : M 22 : M 33 = e 5 : e 4 : 1, M 23 ~ M 22 , (222) 

the resulting effective neutrino mass matrix has the following form 

/ ee 2 ^ \ „2 

M e J ! = ee 2 ^ -cre 2 ^ + 2pe^ rje** • (223) 
V V e<* e 2i *J M3 

where M 3 is the heaviest eigenvalue of M v ^rr which is of the same order as M 33 , 
and a is defined as ere 4 = M 22 /M 3 . The parameters rj, p and a are all of 0(1), and 
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r\ and p are determined using the quark masses. The effective neutrino masses form 
the following pattern, 

m ui : m U2 : m„ 3 = e : e : 1 (224) 

with e ~ 0.1, we have (m 2 2 — m£ )/(m^ — 777 2 2 ) ~ 10~ 2 , which is consistent with 
the LMA solution. Nevertheless, it is not clear whether the predicted solar angle is 
consistent with experiment or not. The prediction for U eV3 is of order 0(e). 

An interesting feature of this model is that, with e ~ 0.1, the scale of Mi is 
about 10 9 GeV . The amount of baryonic asymmetry which is very sensitive to this 
scale is given by 

Y B k sign{a) sin(0 - a) x 10~ 9 . (225) 

The result has the promise to be consistent with the observed value, once the 
parameters K, a, 4> and a are fixed. 



5.2. Matsuda, Fukuyama and Nishiura 

Matsuda, Fukuyama and Nishiur d^^ proposed Hermitian textures with four 
in the context of 5O(10) 



o K,v 


^ 




Cu,v I 


Cum 


D u>v ) 





A d , e e ia ™ 




Bd, e Cd, e e ia 








(226) 



(227) 



For a general matrix of the form 

M = I A B C I , (228) 




one can relate its three eigenvalues, mi, m%, and 7713, to the matrix elements by 

DA 2 = -m 1 m 2 m 3 (229) 
BD — A 2 — C 2 — m\mi + m%m?, + m 3 m\ (230) 
B + D = m 1 +m 2 + m 3 . (231) 

Because of the observed fermion mass hierarchy, 1 777,3 1 3* |ttt-2 | 3> |t77i|, it is a good 
approximation to write B = 7772 and D = 777,3 ~~ m i & n d express the mass matrix in 
terms of its three eigenvalues as 



y/— m 1 m 2 

M ~ I yj — rtl\mi 7772 -y/— 77717773 ] . (232) 

— 777 1 777 3 777 3 + 777! 
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To put this idea to work in the context of 5*0(10), we first note that the set of 
equations given in Eq.(160)-(163) can be re-written as, assuming < 120 > is small 



and contribute to M e and M d only, 

M u = Y w < 10 > +eY~ 126 < 126 > (233) 
M d = aY w < 10 > +Y 126 < 126 > +Y 120 < 120 > 

= aM u - (ae - 1)Y~ 126 < 126 > +Y 12a < 120 > (234) 
rM e = aY 10 < 10 > -3Y~ 126 < 126 > +SY 120 < 120 > 

= aM u - {ate + 3)Y~ 126 < 126 > +SY 120 < 120 > (235) 

r ' M Dirac = ylO < 10 > -'3 e Y 126 < 126 > (236) 

sM v<LL = (3Y 126 < 126 > (237) 

sM^rr = 7 r 126 < 126 > (238) 

where parameters a, 0, 7 are ratios of SM Higgs doublet VEVs from different 
5*0(10) representations. The symmetric (anti-symmetric) matrices Y 10 ' 126 (Y" 120 ) 
can be expressed in terms of the symmetric (anti-symmetric) part of the mass 
matrices M Ut d,e,i/ a s 

(1 - ae)Y 10 < 10 > = (M u ) s - e(M d ) s (239) 

Y 126 < 126 > = ~(M d ) s - K{M e ) s (240) 

Y 120 < 120 > = (M d ) a (241) 



where the subscripts s and a refer to symmetric part and anti-symmetric part, 
respectively. One can then solve for the elements in matrices Y~'s in terms of quark 
masses. As the simple approximations, B — 1112 and D — + mi, work well for 
quark masses and CKM matrix elements, in order to obtain viable neutrino masses 
and mixing angles, a deviation must be made in the charged lepton mass matrix, 
B e = TrifiiX +£) and D e = m T + m e — £m M . With £ ~ 0.01, maximal — v r mixing 
angle and the LMA solution can be accommodated. In this model, the allowed 
region for the leptonic CP violating Dirac phase can be obtained. 



5.3. Bando and Obara 

Bando and Obar J i24 l 12 S I pursue along the line of Matsuda, Fukuyama and 
Nishiurdl^l to analyze mass matrices of the type given in Eq. l|226|l . and have a 
detailed analysis on all possible combinations of contributions from either < 10h > 
or < 126# > for each non-vanishing entry. In other words, all possible ways the 
CG factor (—3) due to < 126# > can appear in the neutrino Dirac mass matrix 

*a L , 

^•dy *0ii *C V (242) 

*c„ 1 
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where * is either 1 or —3 depending upon whether the coupling is due to < IQ H > or 
< 126# >. They found the following texture has best agreement with experiments 

126 











)■ 


a u 


b u 


c u 




Vo 




1 / 



M u = ( 126 10 10 | = ( a u b u c u | . (243) 
10 126, 

In this case, the Dirac neutrino mass matrix is 

/ -3a„ 

M® irac = -3a u b u c u ) m t (244) 
V c u 

and the right-handed neutrino Majorana mass matrix is generated by the coupling 
to 126 Higgs representation, and is of the form 

/0 r 0\ 

M V . RR = r \m R . (245) 

\ooi/ 

The effective neutrino mass matrix is thus given by 

f T • 

^ 2afcVc 2 c(f + 1) I ^ (246) 



Mf/ f 



f2 



V o c(f + 1) d 

where a = — 3a u , b = b u , c = — 3c u , and d = —3. The typical predictions for this 
type of mass matrices are 

sin 2 26*23 ~ 0.98 - 1 (247) 

tan 2 0i2 - 0.29 - 0.46 (248) 

\9 13 1 - 0.037- 0.038 (249) 

\m V3 1 - 0.053 - 0.059 eV (250) 

\m V2 1 ~ 0.003- 0.008 eV (251) 

0.0006 - 0.001 eV. (252) 



in 



6. 50(10) Models with Family Symmetry in 4-dimensions 

A natural framework to accommodate small neutrino masses is a grand unified the- 
ory based on SO(10) in which a right-handed neutrino in each family is predicted 
and the see-saw mechanism can be implemented naturally. Many other models based 
on 50(10) besides those mentioned in Section 5 have been proposed to accommo- 
date the observed fermion masses and mixing angles. Here we concentrate only on 
models which utilize family symmetry. We classify these models according to their 
intermediate symmetry below S'O(IO) breaking scale and the family symmetry. Dif- 
ferent symmetry breaking pattern of 50(10) have different mass relations among 
the quark and lepton sectors, resulting in different way to generate large leptonic 
mixing angles. 
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6.1. Models with Symmetric Mass Textures 

Symmetric mass textures naturally arise if the 5*0(10) is broken down to the SM 
gauge group through the left-right symmetry breaking chain. Due to its symmetric 
nature, this type of models tend to be more predictive compared to models with 
lop-sided/asymmetric mass texture. 



SU(2) Family Symmetry 



6.1.1. Chen and Mahanthappa 

The model proposed by Chen and MahanthappeP^EHKEIhas SU (2) family symme- 
try. Because S'O(IO) breaks down through the left-right symmetry breaking chain, 
symmetric mass matrices arise. The field content of this model is given by, in terms 
of ,50(10) x SU(2) quantum numbers, 

Matter f ields : 1(16,2), 1(16,1) 

Riggs fields : 5(10, 1), 3(T26, 1) 

Flavonf ields : 3(1, 2), 3(1, 3). 

After the symmetry is broken, the following mass matrices are generated 



M„ 



M d , e = 



(10+) e'* 

(l0+)e(l0+)e 
,(l0+)e'(l0 3 +)e (10+) 

<10 5 7)e' ^ 






(lOs)e' (l,-3)(l26 )e | 

(ior)/ 



e' s 
e' (l,-3)pe0 
1. 



Mr 



where 



Mr 



(10+), M D = (lQi) 
r 2 = (l0t)/(l0t), r 4 = (l0+)/(l0+) p = (m~) / (lOf) 



(253) 
(254) 

(255) 
(256) 



eM and e M are the VEV's accompanying the flavon fields. The mass hierarchy 
arises due to the Froggatt-Nielsen mechanism. The right-handed neutrino masses 
are generated due to 126//, thus R-parity is preserved at all energy scales, and no 
additional proton decay modes are allowed, which is to be contrasted to the case 
when 16/r's are implemented to generate right-handed neutrino Majorana masses 
(see the model of Babu, Pati and Wilczek discussed in the next section.) The right- 
handed neutrino mass matrix is given by 



' ^126 2 °^(5i\ 



M VRR = 



/T26l° 



Si 



1262° 
T26 2 ° 



2 
'0 
2 

■v : (126'° 



62 <126,U 




(257) 
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with Mr = ( 12Q 1 V The effective neutrino mass matrix is of the following form 



where t < 1. This model can accommodate the LMA solution, in addition to LOW 
and "Just SO" VO solutions. Its prediction for U ev s is about 0.15. With 11 input 
parameters, this model predicts 22(+9) masses, mixing angles, and CP violating 
phases for quarks and leptons (and right-handed neutrinos). 

We note that in this model, the mass matrices, M ULR , M VRR and M^e//, have 
identical form. In other words, the texture considered in Ea. (|86|) is invariant under 
the se e-saw mechanism. The form invariance also occurs in a model of neutrino 
mixing^^J which uses different texture. 

6.2. Models with Lop-sided/ Asymmetric Mass Textures 

Models having SU(5) as the intermediate symmetry have lop-sided Yukawa matri- 
ces. This is due to the SU(5) relation, = Mj. This opens up the possibility 
of large leptonic mixing angles due to the large mixing angle in the charged lepton 
mixing matrix. 

U(l) family symmetry 
6.2.1. Babu, Pati and Wilczek 

The model proposed by Babu, Pati and WilczelJi^l utilizes the Abelian U(l) as its 
family symmetry. Because it is based on dimension-5 operators to generate fermion 
masses, its field content, in terms of 50(10) quantum numbers, is somewhat simple 

Matter fields : 16i, 16 2 , 16 3 
Biggs fields: 1(10), 1(16 + 16), 45. 

The mass matrices generated are given by 




(258) 




(259) 




(260) 




(261) 
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-3e +T] 

M e = | 3e+rj' -3e + 77 | • m d . (262) 
3e + 7? 1 

The right-handed neutrino Majorana mass matrix is generated by the effective 
operator, ^16^16^16^16^ and is given by 

(x z\ 

M VRR = y • M R (263) 
\zylj 

and the resulting effective neutrino mass matrix is given by 



n cff 



9e 2 (x-z 2 ) 3e y(-3e z+(~3e+a)x) 3e(xy~(3t+a)(x-z 2 ) 

3e ' y(~3t z+(-3e+a)x) -9e 2 y 2 (3e+a)y(3e z-(-3e+a)x) 

3e(xy-(3e+<j)(x-z 2 ) (3e+a)y(3e z-(-3e+<r)x) (3e+a)(-2xy+(3e+cr)(x-z 2 )) _ 

(264) 

The large atmospheric mixing comes from the effective neutrino mixing matrix, 
by choosing the value of parameter y. As y also gives rise to small V c b value, the 
smallness of V c b and the maximality of the atmospheric mixing angle are thus 
related. This model can only accommodate SMA solution for the solar neutrinos. 
The LMA solution can be obtained if an intrinsic LH neutrino Majorana mass term 
arising from a dim- 7 operat or is included, assuming it dominates over the regular 
Type I seesaw A characteristic of this model is the presence of a new 

prominent dim-5 proton decay mode, p — > fi + K°, in addition to the p — > vK + 
mode. This is a consequence of utilizing the 16# to generate neutrino massesP23 
In this model, the #13 angle is predicted to be about 5.5 x 10~ 4 . As all parameters 
are taken to be real, CP is conserved in this model. With 11 parameters, this model 
accommodates 18(+6) masses and mixing angles for quarks and leptons (and RH 
neutrinos), yielding 7 predictions in accord with the data. 

6.2.2. Albright, Babu and Barr 



The model proposed by Albright, Babu and BarJ ^ 4 ! 1 ^ 1 ! 1 ^ 1 ^ has U(l) x Z 4 x Z 4 
as family symmetry. The model has the following particle content, in terms of 
,50(10) representations, 

Matter fields: 16 x , 16 2 , 16 3 , 2(16 8 16), 2(10), 6(1) 

Kiggsfields: 4(10), 2(16® 16), 1(45), 5(1). 

The mass matrices in this model are generated by the Froggatt-Nielsen mechanism, 
and are given by 

fr, \ I r, 5 5' e^ 

M u = e/3 -m„, M d = S a + e/3 | ■ m d 
\0 -e/3 1 / \<5'e^ -e/3 1 
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V 

M VLR = I e I • m u , M e = I 5 -e I • m d . (265) 
\0-el/ XS'e^a + e 1 / 

The right-handed neutrino Majorana mass matrix is generated by the effective 
operator of the type, -p-16il6jl6//16//, and is given by 

(c rf —ber] an \ 
-ben e 2 -e • A R . (266) 
an —e 1 / 

And the effective neutrino mass matrix that accommodates the LMA solution is 
given by 

e/(a-b) \ 

Mt }S ~ | e/(a - b) -e 2 (c 2 - 6 2 )/(a - bf -be/ {a - b) m 2 JX R 
-be/(a-b) 1 / 

-e \ 

-e 2e ml/X R , (choosing a = 1, 6 = c = 2). (267) 
2e 1 / 

This model can also accommodate the SMA, and "Just So" VO solutions. An 
interesting property of this model is that the large mixing angle in atmospheric 
neutrinos is due to the lop-sided structure of M e , in which e ~ 0.1 and a ~ 1, 
giving rise to a large left-handed mixing angle in the (2,3) block of V c .l. The 
matrix M e is related to Md by the SU(5) relation, M e — Mj , thus such a lop-sided 
structure gives rise to a large mixing angle in the right-handed rotation matrix for 
the down-type quarks, Vd,R, making it un-observable. The large solar mixing angle 
is due to the structure in the (1,2) block of 



-e 
-e 



(268) 



Because the large mixing angles in the atmospheric and solar neutrinos are due to 
different mass matrices, the prediction for U eU3 can be made to be extremely small. 
In this model, \U eU3 \ is predicted to be 0.014. With 10 parameters, this model 
accommodates all the 22 masses, mixing angles, and CP violating phases at low 
energies. 



6.2.3. Maekawa 
MaekawJSl 

has proposed a 5*0(10) model combined with an anomalous U(\)a 
symmetry. This anomalous U(\) symmetry is important for achieving DTS via 
Dimopoulos-Wilczek mechanism and for generating fermion mass hierarchy. The 
U(\)a charge assignments project out terms that destabilize the VEV of 45# , and 
thus guarantee that < 45# > is along the B — L direction, which otherwise can only 
be achieved with fine-tuning and introducing additional Higgs multiplets. These 
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charge assignments also give rise to the fermion mass matrices of the following 
form: 

/A 6 A 5 A 3 \ /A 4 A y / 2 A 3 \ 
M u = I A 5 A 4 A 2 J < H u >, M d = A 3 A 5/2 A 2 < H d > . (269) 
\A 3 A 2 1 / \ A AV2 i J 

The charged lepton mass matrix is M e = Mj ■ rj, where rj characterizes the renor- 
malization effects, and the Dirac neutrino mass matrix is given by 

/ A 4 A 3 A \ 

M mrac = A 7/2 A 5/2 yl/i . y2 ^ < R ^ > 

V A 3 A 2 1 / 

The right-handed neutrino Majorana masses are generated by couplings to 16-dim 
Higgs representations 

/A 6 A 5 A 3 \ 

M V . RR = A 5 A 4 A 2 ■ A 9 (271) 
\A 3 A 2 1 / 

and the resulting effective neutrino mass matrix is 

/ A 2 A 3 / 2 A \ 

M eff = A 3/2 A A l/2 . A - 5?? 2 < H 2 > ^ 
V A AV2 1 J 

The LM matrix is then given by 

/ 1 A 1 / 2 A \ 
U LM = A 1 - 72 1 A 1 / 2 . (273) 
V A A 1 / 2 1 / 

The LMA solution can be accommodated in this model. The prediction for U e u 3 is 
about A which is very close to the current bound from experiment. All order 0(1) 
coefficients in the mass matrices are not specified in this model, thus the validity 
of this model is unclear. 

6.2.4. Shafi and Tavartkiladze 

Shafi and Tavartkiladze proposed a SUSY SO(10) model combined with an anoma- 
lous U(1)h symmetryESl By extending matter content, the U(1)h can account for 
the observed mass hierarchy and mixing angles in the charged fermion sector. In 
this model, the three light lepton families are contained in the additional matter 
fields, 10i, rather than the three usual 16^ in which quarks reside. From this point of 
view, the quark-lepton unification is lost. The bi-large neutrino mixing is achieved 
by introducing three additional SO(10) singlet states. Due to the U(\)h charge 
assignment, only two of these singlets interact with the lepton doublets giving a 
3x2 neutrino Dirac mass matrix, and the neutrino RH Majorana mass matrix is 
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approximately 2x2. The LMA solution can be accommodated in this mdcl. 
U(2) family Symmetry 



6.2.5. Blazek, Raby and Tobe 

The model proposed by Blazek, Raby and TobJSHIISZI has U{2) x U{l) n as its 
family symmetry. It has the following field content, in terms of 50(10) x U(2): 

Matter fields : 1(16,2), 1(16,1), 1(1,2), 1(1,1) 

Riggsfields : 1(10,1), 1(45,1) 

Flavon fields: 2(1,2), 1(1,3), 2(1, l A ). 

The mass matrices of this model have the following form 



Kiep (e + n 2 e)p 
Y,,, = | — (e — n 2 e)p ep erT- | A 
erTQ 1 

' K\ep e + K 2 e 
Yd = — (e — n 2 e) ep 


3ki£CJ 

Yy, t, = I (e + 3K 2 e)u 


3ki£ - 
V, = I (e +3n 2 e) 

erT T a 




\erT-jyUj X 



(274) 



(275) 



(276) 



(277) 



where u> = 2a /(2a — 1) and Tf = (Baryon number - Lepton number) of multiplet 
/. The right-handed neutrino Majorana mass matrix is generated by the effective 
operator of the type, -^-16^16^16^ 16h, and is given by 



(KiS k 2 S ' 
K 2 S S 
o^o 



(278) 



This model can accommodate LMA solution, in addition to SMA and "Just So" 
VO solutions. In the case of LMA solution, its prediction for U e u 3 is 0.049. It has 
16 input parameters which yield 22 masses and mixing angles. 



6.2.6. Raby 

Raby ha s p roposed a 5O(10) model combined with SU(2) x U(l) family 
symmetryp2> in which the ansatz proposed by Frampton, Glashow and Yanagi da 53 
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discussed in Sec. 3 naturally arises. The representations utilized in this model, in 
terms of 50(10) x SU(2), are given as follows: 

Matter fields: 1(16,2), 1(16,1), 3(1,1) 
Biggs fields: (10,1), (16,1), (45,1) 

Flavonfields : (1,1), (1,1 A ), (1,2), (1,2), (1,3), 

in addition to several vector-like Froggatt-Nielsen fields. For the purpose of discus- 
sion, we denote the three 50(10) x SU{2) singlets by Ni (i — 1, 2, 3) which play an 
important role in the neutrino sector, the SU (2) (anti-)doublet flavon fields by (<j>) 
<fi, and the SU(2) singlet flavon fields by 9 and Si (i — 1,2). In the charged fermion 
sector, the Yukawa matrices are given by 

/ e'p-e£\ 
Yu = -ep ep -e A (279) 

V ^ e 1 / 

/ e -e£a\ 
Y d = \ -e e -ea A (280) 

V < e I J 

/ -e'oj |e^w\ 
Y VLR = e'uj 3euj \euj A (281) 
\-3e£cr-3eer 1 / 

/ -e 3eA 
Y e = e 3e 3e A. (282) 
\-3eCcr-3ecr 1 / 

The ansatz of Frampton, Glashow and Yanagida is obtained by considering the 
following superpotential in which the three 50(10) x SU(2) singlets Ni mix with 
the SM singlets in the three families, via the following superpotential 

W v - ^-(N^Wa + N 2 <j> a lQ a + N 3 916 3 ) + -(SiJV? + S 2 N%). (283) 
M 2 

The complete neutrino mass terms are thus given by 

1 , X 

vm v v + VVN + -NM N N (284) 

where m v (oc Y VLR ) is the Dirac mass matrix due to the coupling between v and 
V, which is related to the up-quark mass matrix, and V and M/v are the Majorana 
mass matrices due to the coupling between JVj and u[, and the self coupling of Ni 
in Ea. (|2"g3|) . 

/ 1 0\ 

V = ^ 2 cf? , M N = diag{M 1 ,M 2 ,0). (285) 
M V 9/ 
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After integrating out the heavy SM singlet neutrinos, and Ni, the neutrino 
effective mass matrix is obtained, 

M e J s = m u {V T )- 1 M N V~ 1 ml. (286) 

The key observation is that if one defines 

fa0\ /10\ 
D T = m v {V T )- l M N V = a 6 1, with V = 1 , (287) 

\0 67 Voo/ 

the effective neutrino mass matrix can be rewritten in the following form 

D T M N 1 D (288) 

where 

M?£) <289) 

which is the ansatz proposed by Frampton et al. Note that in this model, because 
the mixing in the charged lepton sector is small, the LM matrix is approximately 
the diagonalization matrix of M%*f , and the prediction of a bi-large mixing pattern 
is not affected. Even though the bi-large mixing pattern can naturally arise in this 
model, the connection between the CP violation in neutrino oscillation and the sign 
of the baryogenesis, which exists in the model of Frampton, Glashow and Yanagida, 
is lost, due to additional CP phases and more complicated structure in this model. 

SU(3) Family Symmetry 



6.2.7. Berezhiani and Rossi 

The model proposed by Berezhiani and RossP^has SU (3) as family symmetry with 
some unspecified discrete symmetries imposed. It has the following field content, in 
terms of 50(10) x SU(3): 

Uiggsfields: 2(10,1), 1(16,1), (16,1), 2(45,1), 1(54,1) 
Ylavonfields: 1(1, 6), 3(1,3), 1(1,8). 

The mass matrices of this model are given as follows: 

£>e l « Ae 

Y u = | y c | , Y d = | Ae 1 " Vet De* B | (290) 

\c 

Y VLR = | y c | , Y e = | k 3 Ae ia y ct De^ faB ) . (291) 

k 2 C kiB D 



(yu o o\ 




y c 


Y d = | 


V y t ) 




(VuO o s 




y c 


), n-( 


V \y tl 
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The right-handed neutrino Majorana mass matrix is taken to be diagonal. Because 
both Y VLR and M URR are diagonal, the effective neutrino mass matrix is also diag- 
onal. Thus the leptonic mixing matrix is proportional to V e> L- Due to the lop-sided 
structure in the (2, 3) block of Y e arising from the SU(5) breaking chain, the max- 
imal mixing angle in atmospheric neutrinos is obtained. Nevertheless, the solar 
mixing angle in this model is very small, which is in the range of SMA solution. 
This model has 14 input parameters in the Yukawa sector. The value of sin 2#i3 is 
predicted to be O(10~ 2 ). 



6.2.8. Kitano and Mimura 

Kitano and MimurJl^l propose a 50(10) model combined with SU(3) x U(1)h 
family symmetry. The three families of matter fields form a triplet of SU(3). The 
Yukawa couplings in this model have the form 

i,j=l pl * 

The field $ is a singlet of SU(3), but it has non-vanishing U(V)h charges. Thus its 
VEV XMpi provides the mass hierarchy, if different generations, i, have different 
U(1)h charge, proportional to — Xj. The field £j's are 3 representation of SU(3), 
whose VEV break the SU(3) family symmetry 

< fx >~ f j M„, < £ 2 >~ (l j M*, < & >~ ( 1 j M* (293) 

where M* is the family symmetry breaking scale. With the U(1)h charge assign- 
ment, the mass matrices of the up- and down- type quarks are 

/A 6 A 5 A 3> 

Y u ~ 77I (< 6 >< 6 >< £i >) ( A 5 A 4 A 2 

Va 3 a 2 i 

A 6 A 5 A 3< 

A 5 A 4 A 2 | (294) 
A 3 A 2 1 

Y d = Y e T ~ ^(< 6 >< 6 >< 6 >) 



(295) 

These Yukawa matrices give rise to the observed hierarchical masses 

m u : m c : m t ~ A 6 : A 4 : 1 (296) 
md '■ m s : mb = rn e '■ nip. '■ fn T ~ A 4 : A 2 : 1 (297) 
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and the CKM matrix 

/I A A 3 \ 
Vckm ~ -A 1 A 2 . (298) 
V-A 3 -A 2 1 / 

The role of these £ fields is that it gives rise to the bi-maximal mixing pattern in 
the neutrino sector. With the E/(l) charge assignment, the neutrino mass matrix is 

/A 2 A A 3 \ (<£ >\ /A 2 AAA 
(< 6 >< 6 >< £1 >) A 1 A 2 < £f > ~ All. (299) 

\a 3 a 2 aV V<#>/ V A 1 I J 

To implement the above idea in 5*0(10) is not trivial, because all the three families 
must have the same charge assignment under U{1)h- This is overcome by intro- 
ducing many Froggatt-Nielsen fields and diagrams so that the phenomenologically 
viable mass matrices discussed above are reproduced. The neutrino masses are gen- 
erated by a complicated mechanism in which each SM singlet in the 16-dim matter 
field mixes with a Froggatt-Nielsen field which transforms as singlet under S'O(IO). 
Thus, instead of a 6 x 6 neutrino mass matrix, this model has a 9 x 9 neutrino mass 
matrix. The Dirac neutrino mass matrix is different from that of the up-type quark 
because they are due to different Froggatt-Nielsen diagrams. Large mixing angle in 
the atmospheric sector can be accommodated, and the LMA solution can also be 
accommodated. 



6.2.9. Ross and Velasco-Sevilla 

Ross and Velasco-Sevillc^^ utilize SU(3) as the family symmetry in combination 
with S'O(IO). Based on the formulation given in Sec. 2f^the up- and down-type 
quark mass matrices are of the form 

/ A'e 8 Ae 3 Ae 3 \ 
M = -Ae 3 A"e 2 A"e 2 M 33 . (300) 
^-Ae 3 A"e 2 1 + A"e 2 J 

Note the expansion parameter e is different for up- and down-quark sectors. The 
Dirac neutrino mass matrix is of the form 

/ C(e 8 ) e 3 (z + e(x + y)) e 3 (x + e{x - y))\ 
M? irac = I e 3 (z + e(x + y)) e 2 (aw + eu) e 2 (aw-eu) • M?j£ c (301) 
\e 3 (x + e(x — y)) e 2 (aw — eu) 1 / 

where g" 3 ac is the (33) component of the Dirac neutrino mass matrix, and the 
real parameters z,a,w and complex parameters x,y,u are of 0(1). Note that the 
higher order terms in symmetry breaking parameters have been included in M® lrac 
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given above, as they are crucial for getting near bi-maximal mixing pattern. The 
right-handed neutrino Majorana mass matrix is assumed to be diagonal 

/mi \ 

M v ,rr = m 2 . (302) 
V mj 

The mass eigenstates of the effective neutrino mass matrix are given by 

- |rk ~ Ne ;^ (303) 
V z + r z 

v 2 ~ (304) 

\J z 1 + r 2 

^3 - v a (305) 

where £ = Arg(z) —Arg(r) and 7- = y/2(zu — a u wy) / z . Here z^ a is the heaviest mass 
eigen states given by 

v _ ( z + ea;)(t/ r + z^j - q/(i/ T - v^) 
v/(z + e(x + y)f + (z + e(x - y)) 2 

and v\j is orthogonal to v a . As one can see in ^3, ^ contributes equally as v T , 
thus the mixing angle between and v T is maximal. A large solar mixing angle 
is also obtained if tan 2 9\i ~ |^/^| 2 is in the right parameter space. Solutions for 
these parameters that are consistent with other charged fermion masses and CKM 
matrix elements have been found; it can accommodate LMA, LOW and "Just So" 
VO solutions. The CHOOZ angel in this model is predicted to be about U eU3 ~ 0.07. 



6.3. Comparisons of Models and Other Issues 

6.3.1. Distinguishing Models using CKM Unitarity triangle 

To distinguish various existing models, clearly we need more precise results from 
the experiments. The most sensitive test are the three angles in the CKM unitarity 
triangle, as suggested in Ref. 141. To illustrate this point, in Table [SJ predictions 
for the three angles of the CKM unitarity triangle from various models are given. 
Clearly, the predictions from these models are very different. Many models will 
be ruled out as soon as one can pin down the values for these three angles more 
accurately. 

6.3.2. Distinguish Models Using sin 2 6 '13 

The two classes of models discussed in this section predict very different relations 
between U eU3 and Amg/Am 2 tm , as discussed in Ref. 56, 142, 143, 48. The pre- 
dictions for sin 0i3 of various 5O(10) models are summarized in Table|S] In the 
model of Chen and MahanthappapSI the typical value for sin #13 is very close to 
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the sensitivity of current experiments. In this model, the value of U eV3 = sin #13 is 



Table 5. Predictions for the three angles of the CKM unitarity triangle from different 
models. 





sin 2a 


sin 2/3 


7 


Albright-BarJESl 


-0.2079 


0.6428 


64° 


Blazek-Raby-TobJEZl (LMA solution) 


0.94 


0.39 


47° 


Berezhiani-RossCSHl (Ansatz B) 


0.27 


0.75 


73° 


Chen-MahanthappeP^Zl 


-0.8913 


0.7416 


34.55° 


RabjEU 


0.92 


0.50 


71.7° 


Experiments 


-0.95 ~ 0.33 


0.79 ± 0.12 


34° ~ 82° 



Table 6. Predictions for sin #13 of various models. The upper bound from CHOOZ experiment is 
sin 0i3 < 0.24. First nine models use 5O(10). Last two models are not based on 5O(10). 



Model 


family symmetry 


solar solution 


sin 013 


Albright-BariGSSl 


U(l) 


LMA 


0.014 


Babu-Pati-WilczeJE3 


U(l) 


SMA 


5.5 x 10" 4 


Blazek-Ftaby-TobJSZl 


(7(2) x U(l) n 


LMA 


0.049 


Berezhiani-Rossi 


5(7(3) 


SMA 


o(io- 2 ) 


Chen-MahanthappJ2Hl 


5(7(2) 


LMA 


0.149 


Kitano-MimurJSHI 


5(7(3) x (7(1) 


LMA 


~ A ~ 0.22 


MaekawJiMl 


U(l) 


LMA 


~ A ~ 0.22 


RabjEl 


3x2 seesaw with 5(7(2) jr 


LMA 


~ m„ 2 /2m„ 3 ~ O(0. 1) 


Ross-Velasco-SevilliUni 


5(7(3) 


LMA 


0.07 


Frampton-Glasho-nESl 
-Yanagida 


3x2 seesaw 


LMA 


~ m„ 2 /2m„ 3 ~ O(0. 1) 


Mohapatra-ParidJ21 


RG enhancement 


LMA 


0.08 - 0.10 



-Rajasekeran 
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related to the ratio Atoq / Am^ tm as 

U eV3 ~ (Ami/Am^) 1 / 3 . (307) 

Thus as this ratio increases, the value of the angle #13 increases. As the LMA 
solution is the most favored solution, the angle #13 in this model is predicted to be 
very close to the current sensitivity of experiments. We note that #13 of this order 
of magnitude leads to observable CP violation in neutrino oscillation. 

In the model of Albright and Barrp HHSI the relation between Amg/Am^ im 
and sin 0i3 is quite different: as the ratio Amg/Amj (m increases, the prediction 
for \U el/3 \ decreases. This is due to the fact that the large mixing angles in the 
atmospheric sector and solar sector have different origins. The most favored LMA 
solution thus implies that the value of U eV3 is extremely small; a neutrino factory 
is needed in this case in order to pin down its value. 

6.3.3. b — t unification 

In most models, the prediction for nib at the weak scale tends to be higher than 
the experimental observed value, and thus a threshold correction of the size —(15 — 
20) % is needed to bring down nib. Such a large threshold corrections for nib are 
expected due to loop diagrams of SU (2)-singlet bottom squark, the SU (2)-doublet 
third g enera tion squark, gluinos and charginos, as discussed in Sec. 4. Barr and 
Dorsnei^S] suggested that, instead of these threshold corrections to m& being large 
and negative, nib/m T may indeed be smaller than one at the GUT scale, and the 
deviation from the naive b — r unification in SU(5) is due to the large off-diagonal 
element of the charged lepton mixing matrix which also explain the large mixing 
in atmospheric neutrinos in models with lop-sided mass textures. 

6.3.4. CP Violation 

CP violation arises in these models from the complex phases in the VEV's of the 
scalar fields, and from the complex phases of the Yukawa coupling constants. Thus 
they are free parameters in the models. In the quark sector, the complex phase 
is constrained by the masses and the mixing angles, thus definite predictions for 
the three angles (a,/3, 7) in the CKM unitarity triangle can be obtained. In the 
leptonic sector, on the contrary, it is not possible at this moment to obtain definite 
predictions for the three CP violating phases. The situation will be improved once 
the absolute scales of neutrino masses and the three mixing angles are known to 
much better precision. 

7. SUSY SO (10) in Higher Dimensions 

The idea of orbifold (SUSY) GU T's was first proposed by Kawamura to solve the 
double t-trip let splitting problem and later developed by Altarelli and 
FerugicESl an( J Hall and Nomura^H The size of this type of extra dimensions are 
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small, being inverse of the GUT scale, R ~ 1/Mqut- To see how it works, let us 
consider the case with only one extra dimension, which is compactified on a S 1 /Z2 
orbifold. The circle S 1 has radius R and is defined by y — y + 2-kR. Under Z%, y 
is mapped to —y. Thus the physical region can be taken as < y < irR. Various 
components transform under the Z2 symmetry as follows 

A„(x, y) -> A„(x, -y) = PA^x, y)p- 1 (308) 
A 5 {x,y) -» A 5 (x,-y) = -PA 5 (x,y)P~ 1 (309) 
9(x,y)->$(x,-y)=±P*(x,y) (310) 

where A^ (x, y) and A§(x, y) are the components of the gauge fields along the usual 
4Z3 and the 5t/i-dimension, respectively; ${x,y) is a generic matter or Higgs field. 
The generators transform according to the following transformation rules, 

p T ap-l = T a^ p T ap-l = _ T a _ 

Here T a are generators of the residual symmetry group while T a are the broken 
generators. The 5D bulk field can be decomposed into a infinite tower of KK states 

1 00 



^fo R 
-y) = ^E^ ) ^ ! f- (313) 



=1 

The mode (j) + is even under the Zi symmetry, 

P0+(x, y) = <f>+(x, -y) = +0+(x, y). (314) 

After compactification, it has a zero mode <f>+ . The <f>- is odd under Z2 transfor- 
mation, 

P(f>-(x,y) = (f>^(x,-y) = -<j)-(x,y), (315) 

thus it does not have a zero mode, and its n-th KK mode has a GUT scale mass, 
(2n + \)/ttR. For the broken generators, the corresponding gauge bosons thus ac- 
quire GUT scale masses; only those corresponding to the un-broken generators have 
AD zero modes, which are then identified as the gauge fields of the little group. From 
the conditions given in Eq. (|311(l , one sees that this symmetry breaking mechanism 
only works for non-Abelian symmetry, because the generators of an Abelian symme- 
try always commutes with the parity operator, P. In other words, it is not possible 
to reduce the rank of a group by Abelian orbifolding, c and additional U{1) sym- 
metries survive, along with the SM gauge group, if the GUT symmetry has rank 
larger than 4. Breaking these £7(1) symmetries can be achieved by the usual Higgs 
mechanism. In orbifold GUT models, because the GUT symmetry is broken by the 



e In ver y limited cases of orbifold breaking by outer automorphism, the rank reduction may be 
possibleHMUll 
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orbifold boundary conditions, one does avoid the task of constructing symmetry 
breaking scalar potential which is usually non-trivial. 

One should also note that in higher space-time dimensions, supersymmetry is 
enlarged: In 5-D, N = 1 SUSY has 8 super charges; it corresponds to N = 2 SUSY 
from the AD point of view. By having different orbifold boundary conditions for 
different components in a "super" multiplet, one can reduce supersymmetry to AD 
N = 1, similar to the case of GUT breaking. To break both SUSY and the gauge 
symmetry by orbifolding, a larger discrete Abelian orbifold is thus needed when 
building a realistic model. A N — 2 hypermultiplet can be decomposed into two 
N = 1 chiral multiplets; and a N = 2 vector multiplet can be decomposed into a 
N = 1 vector multiplet (denoted by "U") and a N = 1 chiral multiplet (denoted 
by "£"). (See, for example, pp. 348-351 of Ref. 84.) To break N = 2 SUSY down 
to N = 1, we thus require V to be even under the parity, and S to be odd. 

An immediate question one might ask is that how does the Georgi-Jarlskog 
(—3) factor arise? Recall that in AD GUT models, the GJ factor arises as the 
Clebsch-Gorden coefficients associated with the VEV's of scalar fields along certain 
symmetry breaking directions. In orbifold GUT scenario, because the GUT sym- 
metry is broken by orbifold boundary conditions, the GJ factor must arise in some 
other way. 

7.1. Breaking SUSY 50(10) by Orbifolding 

Several orbifoldings have been found to break 50(10). The number of extra 
dimensions that has been considered is either one or two. Some orbifoldings 
break ,50(10) to only its maximal subgroups; others break ,50(10) fully down to 
SU (3) x SU(2) l'xU (1)y X U(l) . Here we summarize various orbifolding that have 
been constructed. 



7.1.1. SUSY 5O(10) m 5D 

Dermisek and Maf^22l consider 5O(10) in 5D and the extra dimension is compact- 
ified on 5 1 jiZ^ x Z 2 ) orbifold. The parities are chosen to be 

P = 4x5^2x2 (316) 

P = diag[-l, -1, -1, 1, 1) <g> 7 2x2 (317) 
and the orbifold boundary conditions are given by 



/ic . T/++ T/++ T/+~ V~ + V — I" V~ + V 

' V (15,14) (1,3,1) V (1,1,3) (6,2,2) ^(15,1,1) (1,3,1) (1,1,3) ^(6,2,2) 

1(1 . TJ + + ff c ff c 1" 

±U H ■ (1,2,2) (6,1,1) (1,2,2) (6,1,1)' 



(318) 



After compactification, the parity Z2 reduces N = 2 SUSY to N = 1 SUSY in 
AD, and the parity Z 2 reduces 50(10) to Gps- The residual symmetry below the 
compactification scale is the Pati-Salam group 5J7(4) C x SU(2)l X SU(2)r on the 
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4D "hidden" brane at fixed point y = irR/2, which is then broken down to Gsm 
by the usually Higgs mechanism, and the symmetry on the "visible" brane at the 
fixed point y = is 5*0(10). 

Kim and Raby-1^ pursue along this line, and analyze the renormalization group 
evolution in this model: The 5D gauge coupling constant unification scale M* (to 
be contrasted with the 4D unification scale Mqut) is found to be ~ 3 x 10 17 GeV 
and the compacti fication s cale is found to be ~ 1.5 x 10 14 GeV. 

Kyae and ShafP^HSSI consider different parities, 

P = diag(I 3x3 ,I 2x2 ,-I 3x3 , -hy.2) (319) 

P = diag(-I 3X 3,l2x2,hx3,-hx2) (320) 

and the Z 2 x Z 2 charge assignments for the components of the 5*0(10) gauge field 
is 

45 v : V( 8il)0 ) V (tX,0) V (l,l,0) ^(3,2,-5/6) ^(3,2,5/6) V (3,l, -2/3) 2V (3,2,l/6) V {lS,l) ■ ( 321 ) 

With these boundary conditions, they are able to break 50(10) down to SU(3) x 
SU(2) L x U(1) Y x [/(!)'. 



7.1.2. SUSY 50(10) in 6D 

Asaka, Buchmuller and Covi et al EESI consider 50(10) in 6D and the extra 
dimensions are compactified on a T 2 /(Z 2 x Z 2 G x Z 2 S ) orbifold. The idea is 
based on the observation that a simple extension of the SM gauge group to 
5t/(3) x SU(2) L x U(1) Y x U(l)' (which has the same rank as 50(10)) is the com- 
mon symmetry subgroup of the Pati-Salam gauge group, 5C/(4) x SU (2) l x SU {2)r 
{Gps) and the Georgi-Glashow gauge group, SU(5) x U(l), (Ggg)- The first parity 
P breaks supersymmetris down to N — 1 in AD, upon compactification on T 2 jZ 2 . 
The other two parities break the 50(10) gauge symmetry, and can be taken to be 



Pgg = 



\ 



a 2 



<T 2 



(j 2 



\ 



Pps = 



v 2 / 



/-co 



V 



-CTQ 



-cr 



-cro 



(322) 



-<ro/ 



in the vector representation of 50(10). At the fixed point of Z 2 G , SUSY Ggg is 
respected; at the fixed point of Z 2 S , SUSY Gps is respected. The charge assign- 
ments for the gauge fields are chosen such that component fields belonging to the 
symmetric subgroup have positive parity and those belonging to the coset space 
have negative parity. At the intersection of two 5D subspaces of the 6D bulk, in 
which Gps and Ggg are un- broken, respectively, the extended SM gauge group, 
5t/(3) x SU(2)l x [[/(l)] 2 is realized. At this intersection (which is also one of the 
fixed points of the orbifold transformations), the electroweak symmetry and the 
additional U(l) are broken by the usual Higgs mechanism. 
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7.2. Fermion Mass Hierarchy in SUSY 50(10) Models in Higher 
Dimensions 

A few mechanisms have been proposed to solve the fermion mass hierarchy prob- 
lem in orbifold GUTs. Some models address the realistic mass relations; the gauge 
symmetry breaking in this type of models is due to both orbifolding and the usual 
Higgs mechanism from which the GJ factor —3 needed to achieve the realistic mass 
relations arises as the CG coefficient associated with the VEVs of some Higgs mul- 
tiplets. Other models make use of mechanisms that are purely higher dimensional, 
e.g. overlap wave function, to generate the mass hierarchy. No mechanisms have 
been found to generate the CG factor by compactification. 



7.2.1. Hall et al 

Hall et a /lliiZI proposed three 50(10) models in 6£>, and the two extra dimensions 
are compactified on T 2 /Z2, T 2 /Zq, and T 2 /{Z 2 x Z 2 ) tori, respectively, in each 
of these three models. In the first model there is N = 1 SUSY in the bulk, and 
the other two models have A 2 SUSY. These models incorporate a mechanism 
proposed by Hall et al | 149 | 158 | ^ S0 \ Ye ^ e fermion mass hierarchy problem in 
which the correct mass relations are generated by mixing the brane localized matter 
fields with additional matter fields that propagate along the fixed line. The generic 
Yukawa interactions for a branc-confincd field, ip(x), and a bulk field, <&{x,y), are 
given as follows: 

C~f %{A 0( 5(y) ? A 3 + A 1( 5(y)V' 2 <I>(y) + A 2( 5( 2 /)V'<I>(y) 2 + A35(y)$(y) 3 + A 4 <I>(y) 3 }.(323) 
The effective Lagrangian below the compactification scale is 

^ 3 + ^7^(0) + y^fo) + ^?o) + (324) 

where V = M string ■ R is the volume factor. Thus different Yukawa coupling con- 
stants in the AD effective Lagrangian have different volume suppression factors, 
depending upon how many bulk fields are involved in the interactions. Therefore 
by having different matter multiplets locate at different locations in the bulk, the 
mass hierarchy can be generated. 



7.2.2. Albright and Barr 

Albright and B ar J - ^ ^ consider 5O(10) in 5D and the fifth dimension is compactified 
on a S 1 /(Z 2 x Z 2 ) orbifold. All quarks and leptons and most of the Higgs fields 
employed in Ref. 131 are confined to the 5O(10) 3-brane; only the 5O(10) gauge 
fields and a 10- and a 45-dimensional Higgs fields are placed in the 5-D bulk. As a 
consequence, most features of the AD modeff^] a l so exist here; the only exception 
is that because the DTS problem is solved by orbifolding, the Higgs superpotential 
in this case is simpler as the terms needed for DTS are absent. The parity Z 2 breaks 
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N = 2 SUSY (from the 4D point of view) down to N = 1 SUSY; Z 2 breaks 50(10) 
down to the Pati-Salam group. The orbifold boundary conditions under Z 2 x Z 2 
are given by 

a c . t/++ t/++ i/++ l/^ — v — I - v — I - v~+ v — 

■ V (15,14) (1,3,1) ^(1,1,3) V (6,2,2) (15,1,1) ^(1,3,1) ^(1,1,3) ^(6,2,2) 
45jJ : ffjyj) ^(1,3,1) ^(1,1,3) ^(6,2,2) ^(15,1,1) ^(1,3,1) ^(1,1,3) ^(6,2,2) (325) 

in • f/ - ' — ^ f/^ — — fi c ~ "t - 

1U H ■ "(1,2,2) (6,1,1) (1,2,2) "(6,1,1)' 

The VEV < 45h > arising from the complete Higgs superpotential on the visible 
5O(10) brane is along the (B — L) direction, thus the GJ factor of —3 remain in 
this model. Because all the matter fields are confined to the AD brane, the Yukawa 
sector of this 5D model is essentially the same as that given in Ref. 131. 



7.2.3. Kitano and Li 

Kitano and lP^H propose a supersymmetric 50(10) model in 5D: the extra di- 
mension is compactified on S\jZ 2 orbifold, which breaks N = 2 SUSY down to 
N = 1, from the 4D point of view. The gauge symmetry 50(10) can be broken 
either by the orbifold boundary conditions or by the usual Higgs mechanism. All 
three families of matter fields along with the gauge fields propagate in the bulk; a 
45#, a pair of 16 © 16 which are needed to break the rank of the symmetry, and a 
Wh are confined to the visible brane. 

In this model, the fermion mass hierarchy is accommodated utilizing the overlap 
between zero mode profiles along the fifth dimension, as discussed in Sec|21 In a 
5D SUSY theory compactified on 5 1 jZ 2 orbifold, the zero mode wave function of 
a 5-D bulk field with a bulk mass term to is localized exponentially as 

fo(y) ~ e~ m *. (326) 

In the 50(10) symmetric limit, all fields in one family must have the same bulk 
mass term, to^, resulting in unrealistic mass spectrum. When the U(l)x subgroup 
of 50(10) is broken by the Higgs mechanism, the VEV of the scalar field < <fi > 
which triggers this breaking also contributes to the bulk mass terms of the matter 
fields. The resulting bulk mass terms are of the form 

to, -> to, - V2g x Q x <(/)>. (327) 

Because different 51/(5) components of 50(10), 1, 5 and 10, have different U{l)x 
charg realistic mass spectrum can be obtained 

/A 6 A 5 A 3 \ /A 4 A 3 A 3 \ /A 2 A A\ 

Y u ~ A 5 A 4 A 2 , Y d = Yj ~ A 3 A 2 A 2 , ~ A 1 1 . (328) 

\A 3 A 2 I / \ A 1 1 J V A 1 1/ 

The lop-sidcdness of Y e thus gives the maximal atmospheric mixing angle and the 
LMA solution to the solar neutrino problem is accommodated. 
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7.2.4. Shafi and Tavartkiladze 

Shafi a nd Ta vartkiladze considered SO(10) in 5D compactified on a S 1 /(Z2 x Z 2 ) 
orbifoldP^] The flavor structure of this model arises at the fixed point which has 
Gps symmetry, thus the mechanism which generates fermion mass hierarchy and 
mixing angles is purely 4-dimensional. By extending the matter content of their 
model and by imposing U(1)h symmetry, the charged fermion masses and mixing 
angles can be accommodated. Bi-large neutrino mixing pattern is achieved by im- 
posing "flavor democracy" in the neutrino sector, which has been made possible 
due to the extension of matter content. 

8. Conclusion 

SUSY GUT is one of the promising candidates for physics beyond the standard 
model: the hierarchy problem is solved, charge quantization is explained, gauge 
coupling constants unification is achieved; as a consequence, a prediction for the 
weak mixing angle sin 2 W is obtained. It provides a natural framework for small 
neutrino masses to arise, and it has the promise for baryogenesis. We have seen 
in this review how the fermion mass hierarchy can arise from a very contrained 
framework of 50(10); this is achieved by imposing family symmetries. As proton 
decay has not been observed, SUSY GUT's in AD are under siege. This situation 
can be alleviated if the SUSY GUT model is constructed in higher dimensions. The 
presence of extra dimensions also provides new ways to understand fermion mass 
hierarchy. On the experimental side, one hopes that more precise measurements for 
the masses and CKM matrix elements will enable us to distinguish these models, 
thus pointing out the right direction for model building. On the theory side, one 
hopes to obtain an understanding of the complicated symmetry breaking patterns 
and charge assignments that are needed in many of these models, which will then 
shed some light on Physics beyond the Standard Model. 
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